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rigid  body  located  at  point  P  in  which  the  vectors  1,  J,  k 
are  fixed 

cosine  (  ) 

beam  torsional  rigidity,  N-m2 

kinetic  energy  minus  potential  energy,  N-m 

unit  vectors  fixed  in  body  B  that,  for  the  initial  position  of 
B,  lie  along  the  axes  x,y,z  and  for  the  final  position  of 
B,  the  deformed  state,  lie  along  the  axes  x3,  y3,  z3 

unit  vectors  along  the  reference  axes  x,y,z  (fig.  1) 

beam  length,  m 

distributed  applied  loads  in  the  x,y,z  directions,  respec¬ 
tively,  N/m 

distributed  applied  moment  about  the  x3  axis  (i.e.,  an  applied 
pitching  moment),  N-m/m 
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the  transformation  matrix  of  direction  cosines  relating  the 
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deflection  of  point  P  in  the  x  direction,  referred  to  as 
axial  deflection,  m 

axial  deflection  of  point  P  due  to  longitudinal  strain  of  the 
elastic  axis,  a  spatial  quasi-coordinate  defined  as 

u  +  (1/2)  f  (v'2  +  w'2)dx,  m 
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deflection  of  point  P  in  the  y  direction,  referred  to  as 
lead-lag  deflection,  m 
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deflection  of  point  P  in  the  z  direction,  referred  to  as 
flap  deflection,  m 

x,  y,  z  a  set  of  Cartesian  coordinates  with  origin  at  the  beam  root,  the 

x  axis  being  along  the  undeformed  beam,  m  (fig.  1) 

x3,  y3,  z3  a  set  of  Cartesian  coordinates  with  origin  at  point  P  with  x3 
remaining  tangent  to  the  elastic  axis  and  y3  and  z.  along 
principal  axes  for  the  cross  section  in  the  deformed  state,  m 

3  a  rotation  (flap)  of  the  body  B  about  the  unit  vector  -  j;  the 

exact  axis  about  which  the  rotation  occurs  is  dependent  on 
the  sequence  of  angular  rotations,  rad 

3Q  angle  with  the  x  -  y  plane  of  the  beam  in  the  example 

problem,  rad 

6  variational  operator;  also  Dirac  delta  function 

6W  virtual  work  per  unit  length  done  by  nonconservative  applied 

loads,  N-m/m 

«♦  vector  of  virtual  rotations  of  the  reference  frame  x3,  y3,  z3 

fixed  in  the  beam  cross  section,  rad 

£  a  rotation  (lead-lag)  of  the  body  B  about  the  unit  vector  k; 

the  exact  axis  about  which  the  rotation  occurs  is  dependent 
on  the  sequence  of  rotations,  rad 

9  a  rotation  (pitch)  of  the  body  B  about  the  unit  vector  i;  the 

axis  about  which  the  rotation  occurs  is  dependent  on  the 
sequence  of  rotations,  rad 

0  pretwist  angle,  rad 

0j  the  third  angle  in  the  series  lag-flap-pitch;  it  is  the  third 

angle  required  to  rotate  body  B  to  the  final  position  after 
first  lag  then  flap  rotations,  rad 

02  the  third  angle  in  the  series  flap-lag-pitch;  it  is  the  third 

angle  required  to  rotate  body  B  to  the  final  position  after 
first  flap  then  lag  rotations,  rad 

k0  the  k  component  of  curvature  of  the  beam  in  the  example  prob¬ 

lem,  m- 1 

r  rotation  per-unit-length  vector  of  the  reference  frame  x3,  y3, 

z3;  the  i  component  is  the  torsion  (angle  of  twist  per  unit 
length)  and  the  j  and  k  components  are  bending  curvatures. 
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<f>  angle  of  twist  due  to  torsional  shear  strain,  a  spatial  quasi- 
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SUMMARY 


Nonlinear  expressions  are  developed  to  relate  the  orientation  of  the 
deformed-beam  cross  section,  torsion,  local  components  of  bending  curvature, 
angular  velocity,  and  virtual  rotation  to  deformation  variables.  These  expres¬ 
sions  are  developed  in  an  exact  manner  in  terms  of  a  quasi-coordinate  in  the 
space  domain  for  the  torsion  variable.  The  entire  formulation  is  independent 
of  the  sequence  of  the  three  rotations  used  to  describe  the  orientation  of  the 
deformed-beam  cross  section.  For  more  common  cases  in  the  literature  in  which 
one  of  the  three  rotation  angles  is  used  as  the  torsion  variable,  the  result¬ 
ing  equations  depend  on  the  choice  of  the  three  angles.  Differences  in  the 
equations,  however,  are  demonstrated  to  be  in  form  only.  The  present  deformed- 
beam  kinematic  quantities  are  proven  to  be  equivalent  to  those  derived  from 
various  rotation  sequences  by  identifying  appropriate  changes  of  variable 
based  on  fundamental  uniqueness  properties  of  the  deformed-beam  geometry. 

This  development  helps  to  clarify  the  issues  raised  in  the  literature  concern¬ 
ing  the  choice  of  the  angles.  The  torsion  variable  used  herein  is  also  shown 
to  be  mathematically  analogous  to  an  axial  deflection  variable  that  has  been 
commonly  used  in  the  literature.  Both  variables  are  quasi-coordinates  in  the 
space  domain  and  have  been  used  in  derivations  based  on  Hamilton's  principle, 
despite  lack  of  rigorous  justification.  Rigorous  applicability  of  Hamilton's 
principle  to  systems  described  by  a  class  of  quasi-coordinates  that  includes 
these  variables  is  formally  established. 


1.  INTRODUCTION 


1 . 1  Background 

To  adequately  model  a  helicopter  blade,  the  deflections  must  be  treated 
as  moderately  large  and  the  equations  of  motion  will  be  nonlinear.  Such  equa¬ 
tions  are  required  for  properly  analyaing  the  stability  and  forced  response  of 
helicopter  rotor  blades  in  hover  and  in  forward  flight.  To  derive  such  equa¬ 
tions,  it  is  necessary  to  first  specify  the  geometry  of  the  beam  both  in  its 
undeformed  state  and  in  its  deformed  state  at  some  particular  instant  in  time. 
For  typical  beam  theories,  this  involves  expressing  the  position  of  a  generic 
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point  P  on  the  elastic  axis,  the  orientation  of  a  frame  consisting  of  the 
axes  normal  to  and  along  principal  axes  for  the  cross  section  at  P,  and  any 
further  deformations,  such  as  warp  of  the  cross  section  out  of  its  nominal 
plane,  to  adequately  specify  the  location  of  every  material  point  in  the  beam 
(ref.  1).  In  the  discussion  of  kinematics  that  follows  we  will  focus  on  the 
variables  used  to  describe  the  location  of  point  P  on  the  elastic  axis  and 
the  orientation  of  the  reference  axes  attached  to  the  beam  cross  section  at  P, 
assuming  that  no  deformation  of  the  beam  cross  section  occurs.  This  would,  in 
general,  imply  a  total  of  six  variables:  three  deflection  variables  for  the 
location  of  P  and  three  angles1  for  the  orientation  of  the  reference  axes 
with  origin  at  P  with  respect  to  a  suitable  space-fixed  reference  frame. 

For  an  Euler-Bernoulli  beam,  however,  the  cross  section  is  assumed  to  remain 
normal  to  the  beam  elastic  axis  during  deformation.  Thus,  two  of  the  three 
angles  can  be  eliminated  by  expressing  them  in  terms  of  derivatives  of  the 
deflection  variables.  It  is  then  necessary  to  express  all  the  kinematical 
quantities  in  terms  of  the  four  remaining  variables. 

It  is  clear  that  a  material  point  can  only  occupy  one  position  at  a  time. 
The  beam  geometry,  by  which  we  refer  collectively  to  the  position  of  the 
material  points  along  the  elastic  axis  and  the  orientation  of  the  cross  sec¬ 
tions  in  terms  of  direction  cosines  for  the  reference  axes  at  P,  is  then 
unique  at  a  specific  time.  When  the  beam  geometry  is  specified,  so  is  the 
strain,  the  torsion  (or  angle  of  twist  per  unit  length) ,  and  the  components  of 
bending  curvature  along  the  principal  axes  of  the  cross  section  at  P,  since 
these  also  are  geometric  quantities.  Although  the  values  of  these  quantities 
are  unique  at  a  particular  instant  in  time,  the  mathematical  variables  used  in 
defining  the  geometry  are  not  unique.  In  reference  2,  Reissner  shows,  for 
example,  that  it  is  possible  to  formulate  a  large-deflection  beam  theory  com¬ 
pletely  independent  of  the  choice  of  deformation  variables.  These  equations 
are  an  extension  of  the  Kirchhof f-Love  equations  (ref.  1)  and  illustrate  the 
obvious  fact  that  the  geometric  quantities  describing  the  beam  deformation 
under  a  given  load,  although  the  expressions  may  look  different  for  different 
choices  of  variables,  are  determined  by  the  loads  alone.  There  are  many  possi¬ 
bilities,  for  example,  in  the  choice  of  angles  used  to  relate  the  orientation 
of  the  reference  axes  at  P  to  a  space-fixed  Cartesian  frame  (fig.  1).  In 
general,  any  change  in  orientation  can  be  described  by  three  successive  rota¬ 
tions  about  three  distinct,  well-defined  axes.  These  sequential  rotations  can 
either  occur  about  space-fixed  axes,  or  about  axes  that  are  fixed  in  an  updated 
reference  frame  resulting  from  a  previous  rotation.  The  angles  themselves  are 
not  unique,  but  values  of  the  resulting  direction  cosines  that  define  the 
orientation  are  unique. 


JIt  is  possible  to  use  variables  other  than  angles  to  describe  changes  in 
orientation  such  as  Rodriguez  parameters  (ref.  2).  In  this  paper,  however, 
we  will  restrict  the  discussion  to  the  use  of  angles. 
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Figure  1.-  Undeformed,  pretwisted  beam  and  reference  Cartesian  axes  x,  y,  z. 


1.2  Discussion  of  Previous  Work 

With  the  assumption  that  the  beam  cross  section  remains  normal  to  the 
elastic  axis  during  deformation,  the  exact  expressions  describing  the  deformed 
beam  geometry  have  been  developed  and  applied  to  a  derivation  of  the  equations 
of  motion  for  a  helicopter  blade  (refs.  3-5).  In  reference  3,  the  exact  direc¬ 
tion  cosines  for  the  reference  axes  with  origin  at  P  were  derived.  These 
direction  cosines  are  invariant  with  the  choice  of  angles  used  in  describing 
the  orientation.  The  direction  cosines  constitute  elements  of  a  transformation 
matrix  [T]  relating  the  reference  axes  with  origin  at  P  and  a  space-fixed 
set  of  Cartesian  axes,  the  x  axis  of  which  is  along  the  undeformed  beam 
(fig-  1). 

In  reference  A,  the  structural  and  inertial  operators  were  derived, 
through  second-degree  nonlinearity,  for  the  rotating  beam.  The  bending  curva¬ 
tures  were  developed  in  the  appendix  of  reference  A.  These  expressions  are 
exact  and  independent  of  the  choice  of  angles  used  to  describe  the  orientation. 
In  reference  5,  expressions  for  the  aerodynamic  loads  are  developed  for  the 
rotating  beam.  To  derive  expressions  for  the  aerodynamic  loads  one  must  make 
use  of  components  of  each  of  the  following  vector  quantities  along  the  refer¬ 
ence  axes  with  origin  at  P:  the  relative  velocity  of  point  P  with  respect 
to  the  fluid,  the  angular  velocity  of  a  rigid  body  B  in  which  are  fixed  the 
reference  axes  with  origin  at  P,  and  the  virtual  rotation  of  B  (fig.  2). 

The  first  of  these  three  items  depends  only  on  the  [T]  matrix,  but  the  latter 
two  must  be  developed  separately.  Only  approximate  forms  of  these  quantities 
are  given  in  reference  5  and  the  exact  expressions,  developed  herein,  have  not 
previously  been  published. 

Torsion  deformation  in  references  3-5  is  described  by  use  of  a  quasi¬ 
coordinate  in  the  space  domain,  the  elastic  component  of  the  angle  of  twist. 

The  variable  is  referred  to  as  a  quasi-coordinate  because  it  is  related  to 
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angles  through  integrals  that  cannot  be  evaluated  in  closed  form.  Hence,  all 
three  angles,  no  matter  how  they  are  chosen,  can  be  eliminated  from  the  deri¬ 
vation  for  an  Euler-Bernoulli  beam.  Although  the  definition  of  the  torsion 
variable  is  mathematically  explicit  in  references  3-5,  it  is  never  stated 
therein  that  it  is  a  quasi-coordinate.  This  may  account  for  some  of  the  dis¬ 
cussions  that  have  appeared  in  the  literature  (refs.  6-8). 

To  better  understand  these  discussions  it  is  helpful  to  consider  a  simple 
approximation.  Sometimes  an  elastic  beam  is  approximated  by  a  system  of 
spring-restrained  rigid  bodies,  the  simplest  of  which  is  a  single  rigid  body, 
hinged  at  a  point  corresponding  to  the  beam  root  to  provide  three  angular 
degrees  of  freedom.  Now,  for  a  rigid-body  system  such  as  this,  the  arrangement 
of  hinges  presents  the  analyst  with  a  choice  of  infinitely  many  distinct  physi¬ 
cal  spring  arrangements  —  each  of  which  constitutes  a  distinct,  well-defined 
physical  system.  For  example,  in  reference  3  differences  in  stability  between 
two  of  these  systems  are  identified.  They  each  have  only  two  degrees  of  free¬ 
dom:  flap  (out  of  the  plane  of  rotation)  and  lead-lag  (nominally  in  the  plane 
of  rotation).  The  only  difference  between  the  two  systems  is  the  hinge 
arrangement.  By  definition,  lag-flap  means  that  the  flap  hinge  leads  or  lags 
with  the  blade;  flap-lag  means  that  the  lead-lag  hinge  flaps  with  the  blade. 
Thus,  lead-lag  motion  is  precisely  inplane  only  for  the  lag-flap  arrangement; 
hence,  the  word  nominally  is  used  in  the  above  definition.  For  small  angles, 
the  equations  are  linear  and  the  two  systems  have  identical  equations  of 
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motion.  Any  degree  of  nonlinearity  produced  by  moderately  large  angles  will, 
however,  produce  differences  in  the  blade  pitch  orientation  for  the  two  sys¬ 
tems  for  given  flap  and  lead-lag  deflections.  These  differences,  although  not 
large,  produce  different  dynamic  behavior  and  stability  for  the  two  systems. 

It  is  clear,  however,  tliat  an  elastic  beam  with  a  particular  set  of  end  con¬ 
straints  is  only  one  physical  system  and  its  kinematics,  stability,  and 
dynamic  behavior  are  not  subject  to  the  analyst's  choice  of  angles  used  in  the 
derivation  of  equations  that  describe  the  system.  While  this  may  seem  clear, 
there  has  been  considerable  discussion  and  some  misunderstanding  of  this 
issue  in  the  literature  (refs.  6-8). 

The  role  of  hinge  arrangement  in  determining  coupled  flap-lag  stability 
of  a  rotating,  centrally  hinged,  rigid-beam  model  of  a  helicopter  rotor  blade 
in  hover  and  in  forward  flight  is  further  addressed  in  reference  6  where  both 
lag-flap  and  flap-lag  systems  are  studied.  The  arrangement  of  the  hinges  is 
shown  to  influence  the  stability  boundaries  for  both  hover  and  forward  flight. 
In  reference  7  two  different  sequences  of  rotations  are  used  to  describe  the 
orientation  of  the  cross  section  of  the  deformed  elastic  beam  and  second-degree 
nonlinear  expressions  are  developed  for  torsion,  the  bending  curvatures,  and 
the  [T]  matrix.  These  expressions  are  used  in  reference  8  to  derive  two  sets 
of  equations  of  motion  for  a  rotating,  elastic-beam  model  of  a  helicopter  blade 
in  hover  and  in  forward  flight.  In  the  two  sets  of  equations,  the  torsion 
deformation  is  described  by  the  appropriate  angles  of  the  two  sets,  which  dif¬ 
fer  from  each  other  and  from  the  variable  of  references  3-5  as  well.  The 
differences  in  the  variables  are  not  mentioned  in  references  7  and  8,  and  in 
reference  8  the  fact  that  the  formulations  must  be  equivalent  is  only  conceded 
to  be  an  unproven  possibility.  The  following  conclusions  are  drawn  in  refer¬ 
ences  6-8:  (1)  the  transformation  sequence  used  in  deriving  the  equations 

influences  the  stability  and  dynamic  behavior  of  an  elastic  beam  (refs.  6,  8); 
and  (2)  the  work  of  references  3-5  is  incorrect  in  some  respects  (refs.  7,  8). 

Because  of  the  importance  of  the  fundamental  geometric  relationships  of 
slender  beams  to  the  nonlinear  analysis  of  rotor  blades,  and  because  of  the 
apparent  differences  that  exist  between  references  3-5  and  references  6-8,  a 
detailed  exposition  of  this  subject  has  been  undertaken  in  this  report.  The 
influence  of  the  choice  of  angles  is  thoroughly  investigated  in  this  report. 

The  fact  that  any  of  the  hinged-beam  approximations  can  be  analyzed  with  any 
choice  of  angles  will  be  demonstrated.  Although  the  hinge  arrangement  is  cru¬ 
cial  to  accurate  analysis,  the  choice  of  angles  used  in  modeling  is  immaterial 
(ref.  9)  as  long  as  singularities,  discussed  in  section  2,  are  avoided 
(ref.  10).  Since  an  elastic  beam,  unlike  the  rigid  beam  with  its  many  possible 
hinge  arrangements,  is  a  single  physical  system,  we  should  not  expect  the 
choice  of  angles  to  influence  the  outcome  of  the  analysis  at  all.  Otherwise 
a  completely  intrinsic  formulation  like  the  Kirchhof f-Love  equations  (ref.  1) 
or  their  recent  extension  by  Reissner  (ref.  2)  would  be  impossible.  Finally, 
the  work  in  references  3-5,  while  essentially  correct,  does  suffer  from  lack 
of  clarity  in  some  areas  and  it  is  hoped  that  this  report  will  clarify  the 
earlier  work  as  well  as  the  ensuing  discussions  in  references  6-8.  Part  of 
the  difficulty  in  reconciling  the  results  of  other  investigations  with  refer¬ 
ences  3-5  stems  from  the  nature  of  the  torsion  variable  used  therein.  The 
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properties  of  this  variable  need  to  be  clarified,  especially  with  respect  to 
its  use  in  conjunction  with  Hamilton's  principle. 


1.3  Procedure 

In  this  report,  the  appropriate  relations  describing  the  geometry  of  the 
deformed  beam  are  rederived  and  expanded  to  include  the  components  of  angular 
velocity  and  virtual  rotation.  The  appropriate  changes  of  variable  are  estab¬ 
lished  to  prove  that  the  two  formulations  in  references  7  and  8  are  equivalent 
to  each  other  and  equivalent  to  the  approximate  form  of  the  exact  expressions 
derived  herein  and  in  references  3-5.  A  thorough  comparison  with  the  results 
of  other  work  is  undertaken  along  with  a  discussion  of  the  geometric  nature  of 
the  different  torsion  variables  in  references  3-8.  The  consistent  reductions 
of  the  kinematical  equations  for  infinite  torsional  rigidity  and  for  the 
hinged,  rigid-beam  approximations  are  also  examined.  Finally,  the  application 
of  Hamilton's  principle  to  systems  described  by  the  torsion  variable  of  refer¬ 
ences  3-5  is  rigorously  established. 


2.  THE  LARGE-DEFLECTION  GEOMETRY  OF  AN  EULER-BERNOULLI  BEAM 


In  this  section  the  large-deflection  geometry  of  an  Euler-Bernoulli  beam 
is  derived.  For  the  purpose  of  discussing  the  geometry,  we  will  not  consider 
warp.  Without  warp,  cross  sections  remain  plane  and  normal  to  the  deformed- 
beam  elastic  axis  in  an  Euler-Bernoulli  beam.  Thus,  the  geometry  is  completely 
determined  by  three  deflections  u,  v,  w  along  the  x,  y,  and  z  axes  (fig.  2) 
and  some  appropriate  measure  of  the  elastic  twist.  We  must  then  express  all 
other  geometric  quantities  in  terms  of  these  four  deformation  variables.  The 
orientation  of  the  reference  axes  with  origin  at  P  with  respect  to  the 
space-fixed  axes  is  given  by  elements  of  the  [T]  matrix.  Also,  expressions  for 
the  torsion  and  components  of  bending  curvature,  angular  velocity,  and  vir¬ 
tual  rotations  are  developed. 

We  consider  an  initially  straight  beam  segment  with  the  elastic  axis 
along  the  x  axis,  illustrated  in  figure  1.  In  its  undeformed  state  the  beam 
is  pretwisted  so  that  the  principal  axes  of  the  cross  section,  denoted  by  y3 
and  z3,  are  at  an  angle  9t(x)  with  the  y  and  z  axes,  respectively,  and 
0t(O)  is  denoted  by  0Q.  The  cantilever  root  boundary  condition  is  considered 
here  for  the  purpose  of  illustration.  In  principle,  any  set  of  boundary  con¬ 
ditions  could  be  treated  by  slightly  modifying  the  method  that  follows. 

When  the  beam  bends  and  twists  to  some  deformed  configuration,  as  in 
figure  2,  the  geometry  can  be  completely  described  by  the  deflections  u,  v, 
w  of  the  point  P  on  the  elastic  axis  and  by  the  direction  cosines  of  the 
reference  axes  x3,  y3,  and  z3  where  x3  is  tangent  to  the  elastic  axis  of 
the  deformed  beam.  The  matrix  of  direction  cosines  of  x3,  y3,  z3  with 
respect  to  x,  y,  z  constitutes  a  coordinate  transformation  matrix  denoted 
by  [T]  between  the  unit  vectors  i,  j,  k  associated  with  the  local  reference 
axes  x3,  y3,  z3  of  the  deformed  beam  and  the  unit  vectors  I,  J,  K  asso¬ 
ciated  with  the  space-fixed  axes  x,  y,  z.  The  relationship  is  given  by 


properties  of  this  variable  need  to  be  clarified,  especially  with  respect  to 
its  use  in  conjunction  with  Hamilton's  principle. 


1.3  Procedure 

In  this  report,  the  appropriate  relations  describing  the  geometry  of  the 
deformed  beam  are  rederived  and  expanded  to  include  the  components  of  angular 
velocity  and  virtual  rotation.  The  appropriate  changes  of  variable  are  estab¬ 
lished  to  prove  that  the  two  formulations  in  references  7  and  8  are  equivalent 
to  each  other  and  equivalent  t  >  the  approximate  form  of  the  exact  expressions 
derived  herein  and  in  referencjs  3-5.  A  thorough  comparison  with  the  results 
of  other  work  is  undertaken  along  with  a  discussion  of  the  geometric  nature  of 
the  different  torsion  variables  in  references  3-8.  The  consistent  reductions 
of  the  kinematical  equations  for  infinite  torsional  rigidity  and  for  the 
hinged,  rigid-beam  approximations  are  also  examined.  Finally,  the  application 
of  Hamilton's  principle  to  systems  described  by  the  torsion  variable  of  refer¬ 
ences  3-5  is  rigorously  established. 


2.  THE  LARGE-DEFLECTION  GEOMETRY  OF  AN  EULER-BERNOULLI  BEAM 


In  this  section  the  large-deflection  geometry  of  an  Euler-Bernoulli  beam 
is  derived.  For  the  purpose  of  discussing  the  geometry,  we  will  not  consider 
warp.  Without  warp,  cross  sections  remain  plane  and  normal  to  the  deformed- 
beam  elastic  axis  in  an  Euler-Bernoulli  beam.  Thus,  the  geometry  is  completely 
determined  by  three  deflections  u,  v,  w  along  the  x,  y,  and  z  axes  (fig.  2) 
and  some  appropriate  measure  of  the  elastic  twist.  We  must  then  express  all 
other  geometric  quantities  in  terms  of  these  four  deformation  variables.  The 
orientation  of  the  reference  axes  with  origin  at  P  with  respect  to  the 
space-fixed  axes  is  given  by  elements  of  the  [T]  matrix.  Also,  expressions  for 
the  torsion  and  components  of  bending  curvature,  angular  velocity,  and  vir¬ 
tual  rotations  are  developed. 

We  consider  an  initially  straight  beam  segment  with  the  elastic  axis 
along  the  x  axis,  illustrated  in  figure  1.  In  its  undeformed  state  the  beam 
is  pretwisted  so  that  the  principal  axes  of  the  cross  section,  denoted  by  y3 
and  z3,  are  at  an  angle  9t(x)  with  the  y  and  z  axes,  respectively,  and 
6t(0)  is  denoted  by  8Q.  The  cantilever  root  boundary  condition  is  considered 
here  for  the  purpose  of  illustration.  In  principle,  any  set  of  boundary  con¬ 
ditions  could  be  treated  by  slightly  modifying  the  method  that  follows. 

When  the  beam  bends  and  twists  to  some  deformed  configuration,  as  in 
figure  2,  the  geometry  can  be  completely  described  by  the  deflections  u,  v, 
w  of  the  point  P  on  the  elastic  axis  and  by  the  direction  cosines  of  the 
reference  axes  x3,  y3,  and  z3  where  x3  is  tangent  to  the  elastic  axis  of 
the  deformed  beam.  The  matrix  of  direction  cosines  of  x3,  y3,  z3  with 
respect  to  x,  y,  z  constitutes  a  coordinate  transformation  matrix  denoted 
by  [T]  between  the  unit  vectors  i,  J,  k  associated  with  the  local  reference 
axes  x3,  y3,  z3  of  the  deformed  beam  and  the  unit  vectors  I,  J,  K  asso¬ 
ciated  with  the  space-fixed  axes  x,  y,  z.  The  relationship  is  given  by 


When  any  two  of  the  axes  about  which  rotations  occur  have  the  possibility  of 
being  coincident  within  the  primary  range  of  interest  of  the  rotations,  singu¬ 
larities  result.  This  is  the  disadvantage  in  using  what  are  commonly  known  as 
classical  Euler  angles  (ref.  10)  since  the  first  and  third  rotations  are  about 
the  same  axis  when  the  second  angle  is  zero.  The  remedy  is  to  have  the  rota¬ 
tions  about  axes  that  do  not  approach  one  another  for  rotations  in  the  neigh¬ 
borhood  of  zero  since  small  deformations  are  a  special  case  of  interest  for  the 
present  problem.  These  modified  angles,  as  used  in  this  development,  are  also 
sometimes  called  Euler  angles  in  the  literature  (e.g.,  ref.  11).  We  assume 
that  the  unit  vectors  i,  j,  k  are  fixed  in  a  rigid  body  B  and  are  initially 
coincident  with  axes  x,  y,  z.  Rotations  5,  0,  and  6  of  B  occur  about 
k,  - j ,  and  i,  respectively,  but  not  necessarily  in  that  order,  so  that  i,  j, 
k  are  finally  aligned  with  the  x3,  y3,  z3  axes.  Therefore,  the  angular 
displacement  0  rotates  the  beam  cross  section  about  what  is,  for  5  =  8  =  0, 
the  elastic  axis;  the  angular  displacement  B  rotates  the  beam  cross  section 
out  of  what  is,  for  5=0=0,  the  horizontal  (x-y)  plane;  the  angular  dis¬ 
placement  5  rotates  the  beam  cross  section  in  what  is,  for  0  =  0  =  0,  the 
horizontal  plane.  The  first  rotation  always  occurs  about  either  x,  y,  or  z. 
The  orientation  of  the  exact  axis  about  which  the  second  rotation  occurs  is 
dependent  on  the  axis  and  magnitude  of  the  first  rotation.  The  third  rotation 
always  occurs  about  either  x3,  y3,  or  z3.  In  the  text  of  this  report,  only 
two  of  the  six  possible  sequences  of  angles  will  be  considered:  5.,  B3,  0j, 
and  02,  52,  ®2*  The  subscriPts  are  added  to  e^nphasi^e  the  fact  that,  for  a 
given  orientation  of  x3,  y3,  z3  with  respect  to  x,  y,  z,  C2  ^  C x »  82  *  Bj, 
and  @2  03.  Since  it  is  clear  that  the  choice  of  angles  will  not  change  the 

outcome  of  an  analysis,  it  is  natural  to  inquire  at  this  stage  why  it  is  neces¬ 
sary  to  consider  two  different  sequences.  It  is  one  purpose  of  this  report  to 
study  the  outcome  of  expressing  the  deformation  kinematics  in  terms  of  the  two 
sequences  and  to  address  the  issues  listed  in  section  1.2. 

The  axes  of  rotation  for  the  5j,  Bj,  0j  sequence  may  be  described  by  a 
rotation  53  about  the  z  axis  resulting  in  a  new  set  of  axes  xj,  yj,  zj, 
followed  by  a  rotation  8j  about  the  -y|  axis  resulting  in  a  new  set  of 
axes  x£,  y\,  z\,  finally  followed  by  a  rotation  01  about  the  x|  axis 

resulting  in  the  x3,  y3,  z3  axes  (see  fig.  3).  The  axes  of  rotation  for  the 

02,  52,  02  sequence  may  be  described  by  a  rotation  82  about  the  -y  axis 
resulting  in  a  new  set  of  axes  x|,  y|,  z3,  followed  by  a  rotation  52  about 

the  Zj  axis  resulting  in  a  new  set  of  axes  x2,  y2,  z2,  finally  followed  by 

rotation  02  about  the  x|  axis  resulting  in  the  x3,  y3,  z3  axes  (see 
fig.  4).  The  final  set  of  axes  must  be  the  same  regardless  of  the  sequence 
of  rotations.  The  matrix  [T]  may  be  formulated  in  terms  of  any  set  of  three 
angles  that  will  rotate  the  axes  x,  y,  z  to  x3,  y3,  z3-  The  discussion 
here  is  limited  to  the  sequences  5j,  S3,  and  62,  52,  02-  However,  all 

six  possible  sequences  of  the  angles  5,  B>  and  0,  as  defined  above,  are  used 
to  develop  the  geometry  in  the  appendix  of  this  report.  The  conclusions  hold 
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for  any  set  of  angles,  regardless  of  definition,  as  long  as  they  will  describe 
the  orientation  of  axes  x3,  y3,  z3< 

The  geometry  of  the  deformed  state  is  now  developed.  The  position  vector 
describing  any  point  on  the  deformed-beam  elastic  axis  can  be  written  as 
(see  fig.  2) 


r  ■  (x  +  u)I  +  vJ  +  v* 

The  unit  vector  tangent  to  the  elastic  axis  of  the  deformed  beam  is 

=  (x  +  u)+I  +  V+J  +  w4* 


(2) 


(3) 


where  r  is  the  curvilinear  distance  coordinate  along  the  deformed-beam  elas¬ 
tic  axis  and  (  )+  =  3/3r(  ).  Since  the  cross  section  remains  normal  to  the 
elastic  axis  during  deformation 

JS-i-TuI  +  T12J  +  T13K  <*) 


where  T^  is  the  element  on  the  ith  row  and  in  the  jth  column  of  [T ] . 
Thus 


■11 


(x  +  u) 


T1 2  -  v+ 


T!3  =  W+ 


(5) 


no  matter  what  angles  have  been  used  to  express  [T ] .  Since  the  [T]  matrix  is 
orthonormal 


and  a  relationship  between 
and  (6) 


T2  +  T2  +  T2  =1 
11  12  13 


x  and  r  is  easily  obtained  from  equations  (5) 


(x  +  u)+  =  /l  -  v+2  -  w+2 


(6) 


(7) 


or 


r’  =  /(I  +  u')2  +  v'2  +  w'2  (8) 

where  (  )'  ■  3/3x(  ).  The  preceding  equations  must  be  combined  with  explicit 
expressions  for  the  elements  of  [T]  to  completely  specify  the  geometry  of  the 
deformed  beam. 


9 


10 


The  virtual  rotations  are  needed  to  express  the  virtual  work  of  applied  exter¬ 
nal  moments  (such  as  the  aerodynamic  pitching  moment).  The  torsion  (or  angle 
of  twist  per  unit  length)  and  the  bending  curvatures  may  be  deduced  with  the 
use  of  Kirchhoff's  kinetic  analog  (ref.  1)  by  replacing  (*)  with  (  )+  in 
equation  (10) 

k  =  c+k  -  B^(c^  J  -  s^  I)  +  0+i 

=  Kji  +  tCjj  +  (14) 

where  k^,  the  torsion,  is 

<i  =  6^  +  S+sg  (15) 

and  the  components  of  bending  curvature  in  the  j  and  k  directions  are 

kj  = 

+  A  a+  <16> 

*k  =  ?ic61cel  +  B1S01 

We  now  define  an  angle  of  elastic  twist  ()>  so  t’-'at 

(6t  +  *)+  -  <i  (17) 

where  0+=  Q}-x+.  We  note  that  0t  and  <f>  can  be  considered  together  because 
each  occurs  about  the  beam  elastic  axis.  This  angle  of  elastic  twist  is  dis¬ 
cussed  in  references  12-15  and  used  in  references  3-5  as  the  torsion  variable. 
From  equations  (5)  and  (9)  we  find  that 

s  *  w+  (18) 

S1 


Then,  we  differentiate  equation  (5b)  and  make  use  of  equations  (5),  (7),  (9), 
and  (18)  yielding 


(19) 


Substitution  of  equations  (18)  and  (19)  into  equations  (15)  and  (17)  yields  an 
expression  for  0"?  in  terms  of  variables  v,  w,  and  <j>  which  does  not  depend 
on  any  of  the  angles  Cj,  8j»  0^ 


0 


+ 

1 


(0t  +  4>)+ 


+ 


V+W+W++N 

i~^) 


(20) 
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THT 


or 


’• '  8<  +  *  *  f  <v++ + 


dr 


(21) 


if  6^0)  -  6t(0)  *  0Q  and  <J>(0)  *  0. 


The  [T]  matrix  may  now  be  expressed  in  terms  of  v,  w,  and  <j>  with  no 
dependence  on  Cj,  0j,  and  fj.  We  make  use  of  equations  (5),  (7),  and  (9) 
obtaining 


[T]  = 


/l  -  v+2  -  w+2 


-v+cA  -  w"*"s  /l  -  v+2  -  w*2  c.  /l  -  v+2  -  w+2  -  v+w+s 


e. 


AT^7 F? 


A  -  w+2 


s  A  -  w+2 

el 


V+S  -w+c.  A  -  v+2  -  w+2  -s .  A  -  v+2  -  w+2  -  v+w+c. 


sr~-  w+2 


A  -  w+2 


AT-  w+2 


(22) 


where  0j  is  given  by  equation  (21).  This  result  for  [T]  was  first  obtained 
in  reference  3  in  a  different  manner.  The  components  of  the  angular  velocity 
vector,  the  virtual  rotation  vector,  and  the  rotation  per-unit-length  vector 
can  now  be  expressed  in  terms  of  v,  w,  and  <j>.  In  doing  this  we  make  use  of 
equations  (5),  (7),  (9),  (11),  (13),  (15),  (16),  and  (17): 


„  .  w+  /•+  .  v+w+w+  \ 

1  1  -  v+2  _  w+2  \  1  -  w+2  ) 


1  /r 

•  4. 

-vr  c. 


s  A  -  w+2 

0 1 


0)  .  =  - 

J  A" 


=  +  (v+  +  ■  —  -  -1-  -  — 

W+2  V  1  -  W  /  A  -  v+2  -  W 


•  4. 


U)k  =  - 


A" 


+  (v+  + 

w+2  \  1  -  w+2/ 


c  A  -  w+2 
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A  -  v+2  -  w+2 


(23) 
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6*,  -  «e.  + - ^ -  L+  + 

1  1  A  -  v+2  -  w+2  V  1  -  w+2/ 


-c0  6w+ 


3  A  -  W+2 

l+U+*L>+\  ®  1 


«♦  =  1  +  (fiv+  +  V+-+<5^y) - - 

i  A  -  w+2  V  1  "  w+2/  A  - 


v+2  -  w+2 


se  6w+ 


c.  A  -  w+2 
+,,+*,.,+\  ®  1 


Sill, 
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L=+  (sv++^£) — - 

-  w+2  \  1  -  w+2/  /— 


v+2  _  w+2 


<±  *  (8t  +  40+ 

\  1  -  W+2/  A  _ 


-w^c. 


A  -  w+2 


K  =  - 


v+2  -  w+2 


Kk  = 


S01  /  ^  v+w+w++\  Cqi 


A  -  w+2 


r+ 2  -  w+2  / 


(24) 


(25) 


The  variables  0^  and  60j  are  obtained  by  differentiation  and  variation  of 
0 i  in  equation  (21).  Equations  (21)  through  (25)  provide  a  complete  and 
exact  description  of  the  deformed-beam  geometry.  These  equations  are  indepen¬ 
dent  of  ?j,  Bj,  and  0j  since  0j  is  expressed  in  terms  of  v,  w,  and  $. 


To  illustrate  the  uniqueness  of  the  geometry  we  present  the  development 
again  with  a  different  sequence  of  rotations  Bj*  ^2»  02‘  We  use 

equations  (1)  and  (5),  which  do  not  depend  on  the  sequence  of  rotations.  The 
[T]  matrix  can  be  determined  from  the  successive  transformation  relating 
x3,  y3,  z3  to  x,  y,  z  using  figure  4. 


[T]  = 


Cr  Cfi 


CB2C02S?2  "  SS2S02 


W 


Cfi  Cr  Sfl  Cfi  -  C0  Sr  Sft 

9 2  C2  02  02  6 2  ^2  °2 


-s„  c„  +  s .  s_  c n  -s„  c  c„  cn  +  s„  s„  s 


B2  02  02  ?2  02 


0 2  C2  B2  02  62  02  c2 


(26) 


The  angular  velocity  vector  is,  from  figure  4 

-  -  ~KJ  +  -  SqV  +  02A 


(27) 


13 


the  components  of  which  are 


“i  ■  e2  -  e2=5_ 


"i  '  -e2ce2cc2  +  szse2 


“k  ‘  ;2c02  + 


The  virtual  rotation  components  are 

Sip  =  50.  -  s  50- 

1  Cj  ty 


s*j  -  ”ce2cc26e2  +  s82S{2 


5*k  *  ce2St2  +  se2cS2tB2  ) 


The  torsion  is 


4*  .4- 

Ki  =  02  “  62S£, 


(28) 


(29) 


(30) 


and  the  bending  curvatures  are 


'  ~62ce2cc2  +  Ve2 

Kk  ‘  4\  +  e+2\% 


(31) 


Substitution  of  equation  (17)  into  equation  (30)  yields  an  expression  for  0+ 

e*  -  (0t  +  4>>+  +  (32) 


Equations  (5)  and  (26)  lead  to 


s.  =  v+ 


w 


v+w+v++ 


/l  -  v+2  -  w+2  (1  -  v+2)/l  -  v+2  -  w+2  J 


(33) 


Substitution  of  equations  (33)  into  equation  (32)  yields  a  relationship 
between  02  and  v,  w,  and  <p 


A'  >  !  ^ 


(34) 


or 


-et  +  *+  r~  4 

c  J  A  -  v+2  -  w+2  V  1  -  v+2/ 


dr 


(35) 


where  02(O)  =  0 1 (0)  =  0Q  and  <J>(0)  =  0.  Equations  (5),  (26),  and  (28) 
through  (31)  can  be  used  to  express  the  exact  deformed-beam  geometry  in  terms 
of  v,  w,  and  <p.  The  [T]  matrix  is 


[T]  = 
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/l  -  V+2 
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A  -  v+2 


(36) 


The  components  of  the  angular  velocity  vector,  the  virtual  rotation  vector, 
and  the  rotation  per-unit-length  vector  can  now  be  expressed  in  terms  of 
v,  w,  and  <(>.  In  doing  this  we  make  use  of  equations  (5),  (7),  (19),  (26), 
and  (28)  through  (31): 
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s  /l  -  v+2 
02 


A  -  v+2  -  w+2  /l  -  v+2  ) 


(37) 
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Sib.  *  68- - 

1  2  A  -  v+*  -  wt* 


*  rr^) 


6i|< 


<5^, 


j '  -(4”+  +  rr^) 


c  /l  -  v+2 
02 


6v+s, 


/l  -  v+2  -  w+2  /l  -  v+2 


A  -  v+2 


<5v+c. 


-  (i-*  4  ^  — . -  +  - - 

\  1  -  V  /  /l  -  v+2  -  w"*"^  A  -  v+2 


(38) 


<±  =  <et  +  4) 


^-(”+++r^) 

/  ++  v+w4'v++\ 

Kk  =  Vw  +  i  _  v+2/ 


c  /l  -  v+2 
®2 


vT's, 


/l  -  v+2  -  w+2  /l  -  V+2 


s-  A  -  v+2 
62 


V^C- 


(39) 


/l  -  v+2  -  w+2  /I  -  v+2 

where  62  and  66  may  be  obtained  directly  from  equation  (35).  Equations  (36) 
through  (39),  written  in  terms  of  82  defined  in  equation  (35),  appear  to  be 
different  from  equations  (22)  through  (25)  written  in  terms  of  81  defined  in 
equation  (21).  It  is  possible  to  demonstrate  equivalence  of  the  two  sets  of 
equations  by  simply  solving  each  of  equations  (21)  and  (35)  for  $  and 
equating  the  expressions  for  <f>.  The  following  first  integral  is  obtained: 


8?  =  01  +  tan' 


I  ( - Jifet - \ 

\  A  -  v"*"2  -  w+2  / 


(40) 


so  that 


/l  -  v"*"2  -  w"*"2  +  v+w+c 
l _ _ 

/l  -  v+2  A  -  w+2 
A  -  v+2  -  w+2  -  v+w+s 

l _ l 

A  -  v+2  A  -  w+2 


(41) 


These  relations  were  also  obtained  in  reference  15.  Equations  (40)  and  (41) 
can  also  be  derived  by  equating  the  expressions  for  T23  in  equations  (22) 
and  (36).  Now,  by  substitution  of  equations  (41)  into  equations  (36) 


A  * 


through  (39) ,  the  results  can  be  shovm  to  reduce  identically  to  equations  (22) 
through  (25).  Thus,  the  sequence  of  rotations  used  in  the  derivation  is 
proven  to  be  immaterial  and  the  two  sets  of  equations  are  seen  to  be  related  by 
a  simple  change  of  variable.  The  fact  is  further  illustrated  in  the  appendix, 
in  which  all  six  possible  sequences  are  examined  and  shown  to  be  equivalent. 


3.  COMPARISON  WITH  OTHER  WORK 


The  results  of  the  development  in  section  2,  except  for  slight  differences 
in  notation,  are  identical  to  the  results  in  references  3-5.  Several  discus¬ 
sions  of  the  material  in  references  3-5  appear  in  the  literature  and  it  is  the 
purpose  of  this  section  to  examine  these  discussions. 

In  references  6-8,  the  [T]  matrix,  the  torsion,  and  the  components  of 
bending  curvature  are  obtained  for  both  the  Bj,  61  sequence  and  the 

82,  4 2’  ®2  sequence  of  rotations.  The  expressions  therein  are  approximated 
to  second-degree  accuracy  in  the  deflections  u,  v,  w  and  61  or  92.  The 
authors  denote  0j  and  02  by  the  same  symbol  even  though  0j  and  02  are 
different  angles.  After  recognizing  this  slight  inconsistency  in  notation, 
however,  it  is  quite  easy  to  verify  that  the  geometric  quantities  expressed 
therein  are  equivalent  to  the  two  formulations  in  the  preceding  section. 

Indeed,  the  development  in  the  preceding  section  shows  that  the  sequence  of 
rotations  used  in  describing  the  orientation  of  the  local  principal  axes  of 
the  beam  cross  section  is  arbitrary.  Thus  the  sequence  of  transformations 
cannot  influence  the  outcome  of  a  stability  or  response  analysis.  The  equa¬ 
tions  appear  to  be  different  but  a  change  of  variable  can  always  be  found  that 
will  transform  the  equations  derived  based  on  one  sequence  to  that  of  any 
other.  The  application  of  this  fundamental  idea  to  the  rigid-beam  approxima¬ 
tion  of  an  elastic  beam  will  be  discussed  in  section  6. 

Let  us  now  examine  the  question  of  whether  the  sequence  of  transformations 
used  in  an  elastic-beam  analysis  can  affect  the  results  of  the  analysis  as 
discussed  in  reference  6.  This  involves  demonstrating  that  the  two  develop¬ 
ments  in  references  6-8  are  actually  equivalent  to  each  other  and  to  that  of 
references  3-5  and  that  the  [T]  matrix  and  kinematic  quantities  developed  in 
references  3-5,  despite  the  apparent  differences  with  the  results  of  refer¬ 
ences  6-8,  are  nonetheless  correct.  It  is  convenient  for  comparison  with 
references  6-8  to  approximate  the  geometric  quantities  obtained  in  section  2 
to  second  order  in  bending  and  torsion  deflections. 

Let  us  consider  the  case  with  0t  =  0  and  v+,  w+,  0  ,  02  and  i)>  of 
0(e)  so  that  e2  <<  1.  Equation  (8)  shows  that 

r'2  =  (1  +  u')2  +  v'2  +  w'2  =  1  +  2e  (42) 

where  e  is  the  component  of  longitudinal  strain  of  the  elastic  axis  from 
Green's  strain  tensor.  In  engineering  beam  theory,  e  is  neglected  with 
respect  to  unity.  Here  we  assume  e  =  0 (e2)  so  that  r'  *  1  +  0(e2)  and 
u'  ■  0(e2);  thus,  derivatives  with  respect  to  r  may  be  replaced  by 
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Amrv 


. _ :  .a*'  ' 


derivatives  with  respect  to  x.  For  comparison  we  write  only  the  T23  and  T32 
elements  of  the  [T]  matrix,  the  torsion,  and  the  bending  curvatures  of  the 
?!»  9j  formulation: 


T23  =  0x  +  0(e3) 

T32  =  -61  -  v'w’  +  0(e3) 
+  v"w*  +  0(e3) 

=  -w"  +  0JV"  +  0(e3) 
Kk  =  v"  +  0  j w"  +  0(e3) 


(43) 


where  equations  (21),  (22),  and  (25)  have  been  used.  For  the  B2>  c2,  02 
formulation  these  quantities  are 


T23  =  02  -  v’w'  +  0(e3) 
T32  =  _02  +  0(e3) 

Ki  =  °2  ~  v'w"  +  °(e3) 

Kj  -  -w"  +  02v"  +  0(e3) 
<k  -  v"  +  02w"  +  0(e3) 


(44) 


where  equations  (35),  (36),  and  (39)  have  been  used.  These  quantities  are 
identical  to  the  corresponding  expressions  obtained  in  references  7  and  8. 
Equation  (40)  yields,  to  second  order,  the  following  relationship  between 
01  and  02: 

02  =  6i  +  v'w’  +  0(e3)  (45) 


which,  upon  substitution  into  equations  (44),  produces  equations  (43)  identi¬ 
cally.  These  expressions,  which  form  the  basis  for  deriving  the  equations  of 
motion,  are  thus  equivalent.  Equation  (45)  can  be  thought  of  as  a  change  of 
variable.  If  it  is  applied  to  the  equations  of  motion  based  on  the  82,  S2,  02 
sequence  (flap-lag-pitch)  derived  in  reference  8,  the  result  will  be  identi¬ 
cally  the  other  set  of  equations  of  motion  in  reference  8  that  are  based  on 
the  8 1 ,  01  sequence  (lag-flap-pitch). 

Equations  (43)  and  (44)  both  reduce  to  the  same  expressions  when  the 
elastic  angle  of  twist  is  introduced,  from  equations  (20)  and  (21) 


18 


Equations  (47)  relate  <#> ,  0j,  and  02  and  may  also  be  regarded  as  changes  of 
variable  relating  the  angles  0j  and  02  to  the  elastic  twist  angle  (j>.  When 
these  changes  of  variable  are  applied  to  the  development  in  reference  8,  the 
two  sets  of  equations  derived  therein  based  on  8^,  01  and  based  on 

&2>  ^2’  (w^en  specialized  to  hovering  flight)  are  found  to  be  equivalent  to 
the  set  of  equations  found  in  references  3-5  (which  considered  only  the  hover¬ 
ing  flight  case).  Because  of  the  nature  of  tj>,  this  change  of  variable  is 
most  easily  applied  to  the  energy  and  virtual  work  expressions,  however,  not 
directly  to  the  partial  differential  equations,  as  shown  in  section  5.  The 
outcome  of  the  analysis  does  not  depend  on  the  sequence  of  rotations  and  the 
work  of  references  3-5  stands  as  correct. 

Other  authors  (refs.  16  and  17)  have  had  difficulty  reconciling  their 
work  with  references  3-5.  In  reference  16  it  is  stated  that  the  results  of 
references  3-5  are  in  error  because  of  an  incorrect  expression  for  torsion,  a 
comment  similar  to  the  conclusions  of  references  6-8.  The  [T]  matrix  in 
references  3-5  is  stated  to  be  incorrect  in  references  16  and  17  because  the 
pretwist  angle  is  grouped  with  the  torsion  variable  rather  than  being  imposed 
before  the  elastic  deformation.  Furthermore,  it  is  stated  in  reference  16 
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that  the  quantity  0£  must  appear  in  [ T ] -  The  pretwist  0t  can  be  grouped 
with  0j,  02»  and  <p  in  the  present  formulation  because  all  of  these  rotations 
occur  about  the  same  axis.2  The  development  in  section  2  of  this  report  is 
exact  for  an  Euler-Bernoulli  beam,  and  0£  terms  do  not  (and  cannot,  because 
of  dimensional  reasons)  appear  in  [T] .  In  contrast  to  references  16  and  17, 
there  is  no  difference  in  the  final  form  of  the  [T]  matrix  that  results  from 
interchanging  the  order  of  pretwist  and  elastic  deformations  (see  eq.  (17)). 

In  both  references  16  and  17,  the  two  sets  of  axes  that  the  [T]  matrix  relates 
are  not  the  same  two  sets  of  axes  that  the  [T]  matrix  in  this  report  (and  in 
refs.  3-5)  relates.  In  references  16  and  17,  the  [T]  matrix  relates  the  local 
beam  cross-section  principal  axes  after  deformation  to  the  local  beam  cross- 
section  principal  axes  before  deformation  and  not  to  the  x,  y,  z  axes  used 
here.  Thus,  for  0t  ^  0  a  direct  comparison  of  the  [T]  matrix  and  other 
geometric  quantities  in  references  16  and  17  with  those  of  references  3-5  and 
the  present  report  is  not  possible. 

In  reference  12,  two  of  the  present  authors  established  that  the  lateral 
buckling  load  for  uniform,  slender  cantilever  beams  was  identical  for  two  sets 
of  equations  —  one  in  terms  of  01  and  the  other  in  terms  of  <f>.  These 
equations  were  different  in  appearance,  yet  they  were  equivalent  and  yielded 
the  identical  buckling  load.  The  development  in  reference  12  in  terms  of  <p 
proceeds  from  the  Kirchhof f-Love  equations  (ref.  1)  which  in  no  way  depend  on 
the  sequence  of  rotations  used  to  describe  the  cross-section  orientation. 
Careful  examination  of  the  different  approaches  noted  in  reference  12  illus¬ 
trates  that  it  is  important  to  understand  that  the  two  geometric  angles, 
denoted  by  0j  and  02  in  this  report,  and  the  elastic  angle  of  twist  ij>  are, 
in  fact,  different  quantities.  Some  of  the  confusion  about  the  influence  of 
different  sequences  of  rotations  evident  in  references  6-8  apparently  stems 
from  their  use  of  the  same  symbol  for  the  three  different  quantities  0j,  02, 
and  <(>. 


In  reference  18,  nonlinear  equations  of  motion  for  an  Euler-Bernoulli 
beam  are  derived  for  the  purpose  of  investigating  stability  of  the  beam  with 
applied  loading.  Classical  Euler  angles  are  used  to  describe  the  cross- 
section  orientation.  The  three  axes  about  which  rotations  occur  for  infini¬ 
tesimal  angles  are  not  mutually  orthogonal  for  classical  Euler  angles,  for 
which  two  rotations  occur  about  the  axis  that  lies  along  what  is  nominally  the 
beam  elastic  axis.  Thus,  the  torsion  variable  t  is  defined  in  reference  18 
as  the  sum  of  the  first  and  third  Euler  angles.  The  torsion  and  the  I 

component  of  the  virtual  rotation  vector  •  I  are  given  to  second 

order  in  reference  18  as 


=  x’  +  (v"w'  -  v'w")  +  0(e3) 

6^^  =  6t  -  (6v'w*  -  v'6w')  +  0(e3) 


(49) 


2This  is  not  true  of  the  other  four  sequences 
however.  There,  angular  displacements  03,  0^,  05 
beam  cross  section  about  the  deformed-beam  elastic 


treated  in  the  appendix, 
and  0g  do  not  rotate  the 
axis. 
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These  quantities,  in  terms  of  0,  and  0,,  are,  from  equations  (12),  (25), 
and  (40) 

«i  =  +  v'V  +  0(e3)  =  0^  -  v'w"  +  0(e3)  \ 

}  (50) 

iSi^j  =  60j  +  v'5w’  +  0(e3)  =  602  -  w'6v'  +  0(c3)  f 

Equations  (49)  and  (50)  can  be  shown  to  be  equivalent  to  second  order  by  use 
of  the  following  changes  of  variable: 

t  -  0j  +  y~~  +  0 ( c 3 )  =  02  -  +  0(e3)  (51) 

Thus,  the  equations  of  motion  derived  with  the  aid  of  classical  Euler  angles 
in  reference  18  can  be  shown  to  be  equivalent  to  those  derived  with  the  rota¬ 
tions  defined  in  this  report.  To  formulate  the  complete,  exact  geometry  using 
classical  Euler  angles  and  the  elastic  angle  of  twist,  is  quite  difficult, 
however,  and  may  not  even  be  possible  because  singularities  appear  in  the 
integrals  when  deformations  are  small  instead  of  large  as  in  the  present  for¬ 
mulation.  Also,  the  angles  used  in  this  report,  unlike  classical  Euler  angles, 
bear  a  strong  resemblance  to  the  geometric  change  in  orientation  during 
moderate  deformations. 

The  results  of  this  section  show  that  the  expressions  for  deformed-beam 
geometry  developed  in  section  2  and  in  references  3-5  are  equivalent  to  a 
variety  of  formulations  found  in  the  literature.  The  results  of  references  3-5 
are  also  confirmed.  The  nature  of  the  torsion  variable  <j>  and  how  a  formula¬ 
tion  in  terms  of  <J>  differs  from  one  in  terms  of  or  02  has  not  been 
treated  in  the  literature  and  will  be  addressed  below. 


4.  NATURE  OF  THE  TORSION  VARIABLES  0j,  02,  AND  <J> 


In  this  section  we  will  examine  the  nature  of  the  torsion  variables 
0j,  02,  and  <(>.  The  kinematical  relationships  that  define  these  variables  are 
used  to  establish  differences  among  them.  An  example  problem  is  introduced  to 
aid  in  illustrating  these  differences.  The  geometric  relationship  between 
0j  and  02  is  examined  by  projecting  the  local  cross-section  principal  axes  on 
the  yz  plane.  Finally,  the  elastic  twist  tp  is  shown  to  be  a  quasi¬ 
coordinate  (ref.  10)  in  the  space  domain  because  it  is  related  to  physical 
angles  through  integrals,  and  its  mathematical  analogy  with  the  axial  deflec¬ 
tion  of  a  beam  due  to  longitudinal  strain  is  established. 


4.1  Kinematical  Differences 

The  elastic  angle  of  twist  $  is  defined  to  be  the  elastic  component  of 
the  integral  of  the  torsion  Thus,  from  equation  (17) 
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♦  -  f  <±  dr  -  9t  (52) 

*0 

The  elastic  torsion  moment  is  then  simply  GJ<{>+  (ref.  4). 

The  angle  $  may  be  thought  of  as  the  beam  twisting  caused  by  a  distrib¬ 
uted,  externally  applied  torsion  moment  acting  about  the  local  x3  axis.  The 
angle  <j>  does  not  completely  define  the  orientation  of  the  reference  axes  at 
point  P  even  when  and  gj  (or  v+  and  w+)  are  known  at  P.  In  addition,  0j 
(or  02)  must  be  calculated  to  obtain  the  complete  description  of  the  geometry 
of  a  deformed  beam.  This  requires  knowledge  of  and  Bj  for  all  x  inboard 
of  P.  Both  0 j  and  02  occur  as  rotations  about  x3,  the  deformed-beam  elas¬ 
tic  axis  (i.e.,  the  i  unit  vector).  By  definition,  0j  and  02  are  the  third 
angles  in  the  two  different  sequences  Cj,  Bj,  01,  and  B2»  C2,  e2 »  resPec~ 
tively,  as  described  in  section  2.  Since  the  final  orientation  is  fixed,  the 
angles  01  and  02  must  differ.  The  fact  that  the  elastic  twist  is  not,  in 
general,  equal  to  either  of  the  angles  0j  and  02  is  simply  because  of  the 
kinematic  effects  of  combined  bending  of  the  beam  in  two  directions,  v  and  w. 

If  the  bending  deflections  and  elastic  twist  angle  are  known,  0j  and  02  can 
be  written  according  to  equations  (21)  and  (35). 

Jr 

-  w*  +  ^£)dr 

/l  .  v+2  _  w+2  V  1  -  W+V 

(53) 

Jr 

-  v+ . (w-+  + 

A  -  v«  -  V  1  -  v+  / 

U 

They  can  be  approximated  quite  well  by  the  simpler  expressions: 

/x 

v'Vdx  +  0(e3) 

u 

(54) 

0 x 

02  =  +  d>  +  J  v'w"dx  +  0(e3) 

1  '0 

Thus,  rotations  about  the  beam  elastic  axis  (geometrically  like  elastic  tor¬ 
sion  rotations)  may  be  produced  by  combined  bending  deflections  alone,  even 
without  elastic  twist  <j>  or  pretwist  0t> 

4.2  Example  Problem 

The  following  example  problem  is  included  to  further  confirm  and  illus¬ 
trate  the  differences  between  0.  and  02  as  well  as  to  demonstrate  the  use  of 
the  geometric  relations  developed  in  section  2.  Consider  a  beam  that  is  first 
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PROJECTION  OF  I 


Figure  6.-  View  of  deformed-beam  elastic  axis  in  the  plane  determined  by 

i  and  j  (view  AA  from  fig.  5). 


The  matrix  of  direction  cosines  at  any  location  is  quite  easily  obtained  from 
figures  5  and  6 


“cos( 
<j>*  sin( 


KQr)cos 


Knr)cos  8 


0 


-sin(<Qr) 

cos(KQr) 


cos(KQr)sin  6^ 
sin(ic0r)sin  8Q 


/I 

J 


(57 


We  wish  to  show  that  given  u,  v,  w  and  <p  for  a  beam,  the  equations  of 
section  2  will  provide  an  exact  description  for  the  deformed-beam  geometry. 
Since  «£  is  zero  then  $  *  0.  The  geometric  angle  0j  can  be  calculated 
directly  from  substitution  of  equations  (56)  into  equation  (21)  yielding 


0j  =  tan-1  [sin(ic0r)tan  0  )  (58) 

Substitution  of  equations  (58)  and  (56)  into  equation  (22)  for  [T]  yields 


[T]  = 


cos(K0r)cos  60 


-sin(ic0r)  cos(KQr)sin  8Q" 


sin(«0r)cos  80  cos(ic0r)  sln(K.r)sin  6 


-sin  g. 


cos  3 


o 


(59) 


which  is  identical  to  the  matrix  of  direction  cosines  defining  the  example 
deformed  beam  in  equation  (57).  It  can  be  verified  by  substitution  of  equa¬ 
tions  (56)  and  (58)  into  the  bending  curvature  formulas,  equations  (25),  that 
the  exact  relations  defining  the  problem  in  equations  (55)  are  recovered. 

Thus  the  formulas  of  section  2  provide  means  of  exactly  determining  the  orien¬ 
tation  of  the  deformed-beam  cross  section  when  the  elastic  bending  and  elastic 
torsion  deflections  are  specified. 


The  same  results  for  the  transformation  matrix  and  bending  curvatures  are 
obtained  if  one  uses  the  expressions  derived  from  the  62,  £2,  02  rotation 
sequence,  equations  (35),  (36),  and  (39).  Proceeding  as  above  but  starting 
with  equation  (35)  we  find 

02  *  0  (60) 

and  from  equation  (36) 


[T] 


cos(KQr)cos  30 
sin(ic0r)cos  30 
-sin  80 


-sin(<0r)  cos(icQr)sin  0“ 

cos(<Qr)  sin(KQr)sin  80 

0  cos  B0 


(61) 


The  transformation  matrix  is  the  same  as  obtained  in  equation  (59)  and  thus 
again  duplicates  the  example  problem  transformation  matrix  given  in  equa¬ 
tion  (57).  The  fact  that  0j  ^  ©2  should  clearly  indicate  that  the  differ¬ 
ences  between  0j  and  02  cannot  be  ignored.  If  ©1  and  02  were  treated  as 
the  same  angle  it  would  not  be  possible  to  correctly  describe  the  geometry  of 
the  example  problem  in  terms  of  both  0j  and  02. 


4.3  Geometric  Relationship  Between  0j  and  02 

It  is  helpful  to  visualize  0j  and  02  for  the  general  case  of  coupled 
bending  deflections.  This  may  be  accomplished  by  considering  the  local 
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Figure  7.-  Projection  of  the  local 
cross-section  principal  axes  on  the 
y-z  plane. 


Figure  8.-  The  projection  of  unit  vec¬ 
tors  j  and  k  on  the  y-z  plane. 


cross-section  principal  axes  pro¬ 
jected  on  the  y-z  plane.  In  fig¬ 
ure  7  the  projection  of  the  y3  and 
z3  axes  (the  unit  vectors  j  and  k) 
are  shown  in  the  y-z  plane  (the 
plane  containing  unit  vectors  J  and 
K).  The  beam  cross  section  appears 
skewed  because  it  is  not  in  the  plane 
of  the  J  and  K  unit  vectors.  The 
projections  of  j  and  k,  denoted  by 
jp  and  kp,  are  given  by 


Ip  =  T22J  +  T23K 
kp  “  T32J  +  T33^ 

By  inspection  from  figure  8 
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sin  a,  =  - 

1  M 
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✓  ‘22 
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+  T 
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sin  a. 


32 


fi 


+  T  i 
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32  '  *33 

Expressions  for  the  matrix  components 
T23  and  T32  in  terms  of  0j  or  02 
can  be  found  in  equations  (22) 
or  (36),  respectively 


t23  .  Se  A  -  »+! 


T32  ‘  -’’o/1  - 


(64) 


Combining  equations  (63)  and  (64)  and  dropping  higher-order  terms,  we  find 
that 


oij  -  6j  +  0(e3)  | 

a2  *  02  +  0(e3) j 


(65) 


and ,  according  to  equation  (45) 

9,-0 


1 


v’w’  +  0(e3) 


(66) 
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These  results  show  that  for  small  angles,  8j  and  02  are  essentially  the  kine¬ 
matic  rotations  of  the  j  and  k  axes,  respectively,  about  the  i  axis.  The 
fact  that  they  are  not,  in  general,  equal  is  because  the  j-k  plane  does  not 
remain  parallel  to  the  J-K  plane. 

Returning  to  figure  5,  the  view  of  the  beam  in  the  direction  of  the  I 
axis  provides  a  clear  illustration  of  the  kinematic  rotation  of  the  beam.  The 
orientation  of  the  projections  of  the  j  and  k  axes  in  figure  5  is  consistent 
with  the  geometric  interpretation  of  61  and  02  in  equations  (65)  and  (66)  and 
with  the  example  calculations  showing  that  0j  =  -v'w’  and  02  =  0.  These 
kinematic  rotations,  devoid  of  elastic  twist,  were  previously  identified  in 
reference  3  where  their  effect  on  the  blade  pitch  angle  was  noted.  Hence,  the 
term  "kinematic  pitch  rotation"  was  used  to  describe  these  effects. 


4.4  Mathematical  Analogy  of  the  Quasi-Coordinates  $  and  ue 


It  is  well  known  that  even  in  the  absence  of  longitudinal  strain  of  the 
elastic  axis,  there  may  be  significant  nonzero  axial  deflection  of  a  beam  in 
bending.  The  above  calculations  show  that  even  in  the  absence  of  torsional 
(shear)  strain,  there  may  be  some  nonzero  pitch  rotation  of  the  cross  section 
of  a  bent  beam.  To  examine  this  analogy  in  detail  we  must  first  develop  the 
appropriate  relations  concerning  axial  deflection.  The  longitudinal  strain  of 
the  elastic  axis  may  be  explicitly  determined  from  the  relationship  of  x  and 
r  given  in  equation  (42) .  Thus 


e  =  u 


(67) 


We  now  let  u  -  ue  +  uj,  where  ujj  is  a  kinematic  axial  deflection  due  to 
bending  only  and  ue  is  the  axial  deflection  due  to  longitudinal  strain  of 
the  elastic  axis.  The  strain  becomes 


(u4  +  ub)2  v'2 


e  =  u;  +  u^  + 


V  w 

-f  - —  +  — 
2  2 


i  2 


=  0. 

Thus 

v'2 

+ 

2  + 

2 

w 


*2 


=  0 


(68) 


(69) 


or 


u£  =  A  -  v’2  -  w’2  -  1  (70) 

The  kinematic  axial  deflection,  sometimes  called  foreshortening  (refs.  19-21), 
is  then  given  by 
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(71) 


=  -  j  (v'2  +  w'2)  +  0(e4) 


Thus 


ue  =  u  +  x  -  J  /l  -  v'2  -  w'2  dx 

"o 

1  Px 

=  u  +I  r  (v'2  +  w'2)dx  +  0(e4)  (72) 

•'o 

Comparison  of  equations  (72)  and  (46)  reveals  that  ue  and  <p  are  mathematical 
analogs.  The  deflection  ue  is  devoid  of  kinematic  shortening  due  to  bending. 
The  elastic  twist  ij>  is  devoid  of  kinematic  pitch  rotation  due  to  bending.  In 
the  time  domain,  quasi-coordinates  (ref.  10)  are  usually  based  on  integrals  of 
either  components  of  velocity  or  angular  velocity.  By  Kirchhoff's  kinetic 
analog  (ref.  1)  the  angular  velocity  components  carry  over  to  bending  curva¬ 
ture  and  torsion  relations  while  the  velocity  component  along  the  elastic  axis 
can  be  shown  to  correspond  to  a  longitudinal  strain  approximately  equal  to  e. 
To  obtain  the  latter  correspondence,  the  velocity  of  point  P  must  be  written 
subject  to  the  constraint  that  the  velocity  vector  is  tangent  to  the  beam 
elastic  axis.  The  component  of  strain  from  the  analog  is  the  difference  of 
the  velocities,  for  (’)  *  (  )+,  with  and  without  elongation  of  the  elastic 
axis.  Hence,  both  ue  and  <f>  are  denoted  as  quasi-coordinates  in  the  space 
domain  because  they  are  related  to  the  physical  displacement  u  and  angles 
and  02  through  integrals  and  because  their  derivatives  have  velocity  and 
angular  velocity  components  as  kinetic  analogs. 


5.  QUASI-COORDINATES  IN  THE  SPACE  DOMAIN  AND 
HAMILTON’S  PRINCIPLE 

The  angle  of  elastic  twist  tf>  has  been  shown  in  section  4  to  be  a  quasi¬ 
coordinate  in  the  space  domain  similar  to  the  elastic  component  of  axial 
deflection  ue.  Although  several  derivations  appear  in  the  literature,  making 
use  of  at  least  one  of  these  variables  (refs.  4-6,  8  and  19-21),  it  is  well 
known  that  the  standard  form  of  Lagrange's  equations  (and  thus  Hamilton's 
principle)  is  not  adequate  when  quasi-coordinates  in  the  time  domain  are  used 
(e.g. ,  see  ref.  10).  It  is  natural  to  ask  if  there  are  similar  obstacles  when 
the  above  quasi-coordinates  are  used  in  the  standard  form  of  Hamilton's 
principle.  To  answer  this  question  we  will  consider  both  a  general  development 
and  a  simple  example  derivation. 
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S.l  General  Development 

We  will  first  write  symbolically  the  standard  extended  form  of  Hamilton's 
principle  for  nonconservative  systems  using  u,  v,  w,  and  0j.  For  simplicity, 
we  will  consider  only  a  linearized  stability  analysis  so  that  time  dependence 
is  factored  out  of  every  term.  We  then  will  rigorously  take  the  variation. 

The  change  of  variable  will  be  made  to  variables  ue,  v,  w,  and  $  after  the 
variation.  The  resulting  expression,  if  there  are  no  obstacles,  should  corre¬ 
spond  identically  to  the  expression  obtained  by  starting  the  formulation  in 
terms  of  ue,  v,  w,  and  4>  and  then  taking  the  variation. 

The  standard  extended  form  of  Hamilton's  principle  for  an  unpretwisted, 
nonwarping  Euler-Bernoulli  beam  with  variables  u,  v,  w,  and  ©1  is  given  by 

5  J  J  H^u’  u’’  “»  v>  v'»  v"'  w*  w'»  w">  ei»  ei*  dt 

*t:  J0 

*2  a  L 

J  6W  dx  dt  =  0  (73) 


where  H  is  the  kinetic  energy  minus  the  potential  energy  and  6W  is  the 
virtual  work  of  the  nonconservative  aerodynamic  forces.  To  simplify  the  fol¬ 
lowing  development,  equation  (73)  is  written  for  a  linearized  stability  analy¬ 
sis,  with  time  dependence  factored  out  of  every  term  or  a  static  analysis 
without  time  dependence.  The  variation  of  H  when  expanded  yields 
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601 


0  (74) 


where  the  virtual  work  has  been  expressed  in  terms  of  the  generalized  applied 
loads  and  =  60  j  +  w'5v'  to  second  order.  The  forces  are  resolved  along 

the  x,  y,  z  axes,  respectively,  and  the  pitching  moment  acts  about  the  x3 
axis.  We  now  make  the  change  of  variable  that  transforms  equation  (74)  to  a 
different  expression  which  involves  ug  and  The  relationships  between 
u,  0j,  and  ug,  $  are,  to  second  order,  from  equations  (46)  and  (72) 
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Differentiation  of  equations  (75)  yields 


,  ,  v1 2  w'2 

u  *  ue  -  2  "  2 


9 1  -  <f>’  -  v’V 


(76) 


We  now  consider  H  «  H(u  ,  u',  v,  v',  v",  w,  w’,  w",  <j>,  <|> * )  and  note  that 


e  ’  e 
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It  should  be  noted  that  the  above  variation  in  terms  of  4>  (or  ue)  could  not 
have  been  made  if  d> *  (or  Ug)  were  also  an  explicit  function  of  9j  (or  u) , 
as  might  occur  for  general  quasi-coordinates.  Substitution  of  equation  (76) 
into  equation  (74)  yields 
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+  w  J  v  6w  dx  +  ^viw  +  5w  +  60J  +  66.,  +  Lx6u  +  Ly6v 

•'n 


j  +  w’6v') 


dx  =  0  (78) 


+  LzSw  +  *1^(66 

Up  to  now  the  change  of  variables  in  Hamilton's  principle  is  equivalent  to  a 
direct  substitution  for  ue  and  (j)  into  the  partial  differential  equations 
corresponding  to  equation  (74).  The  single  underlined  terms  appear  in  the 
v  and  w  equations  from  ue,  and  the  double  underlined  terms  are  from  4>.  Not 
all  of  these  terms  are  necessarily  visible  at  this  step,  however,  depending  on 
the  order  of  terms  that  are  retained.  The  example  given  below  should  clarify 
the  nature  of  these  terms.  We  now  substitute  for  66,,  68j,  6u,  and  6u'  in 
equation  (78) 

66  j  =  6<J>  -  J  (w'6v"  +  v"6w')dx 

Jo 

68 J  =  6<}>*  -  w'6v"  -  v"6w' 

6u  =  5ug  -  /  (v'6v'  +  w'6w')dx 

*0 

5u'  =  6ug  -  v’6v'  -  w’ 6w' 

All  of  the  underlined  terms  cancel  out  leaving 
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(80) 
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This  is  precisely  the  same  expression  that  would  have  resulted  from  starting 
with  the  quasi-coordinates  and  taking  the  variation  directly.  Since  the  vir¬ 
tual  work  term  involves  only  a  straightforward  change  of  variable,  Hamilton's 
principle  is  valid  when  written  in  terms  of  quasi-coordinates  ue  and  <(>.  Thus 
the  application  of  Hamilton's  principle  with  variables  ue  and  <f>  reduces  to 


6  J  H  dx  +  f  (Lx6ue  +  Ly<5v  +  L-,<5w  +  Mv  6<p)dx 
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dx  (81) 


'X  ^  v  \  'X 

where  the  following  identity  was  used: 

Lx  L  L 

f  f (x)  C  gCx^dXj  dx  =  f  g(x)  f  f (Xj)dXj  dx 

J0  • 0  Jo 


(82) 


It  is  observed  that  the  use  of  the  quasi-coordinates  ue  and  <)>  has  the  effect 
of  replacing  the  terms  underlined  in  equation  (78)  with  integrals  of  the 
applied  loads  from  both  the  u  and  equations. 


5.2  Example  Derivation 

An  actual  derivation  of  equations  will  be  helpful  in  understanding  the 
results  from  the  general  development  above.  Consider  a  rotating  cantilever 
beam  that  undergoes  coupled  axial  deflection  and  flap  (out-of-plane)  bending 
only.  For  simplicity  we  consider  only  the  static  equilibrium  equations  (no 
time  dependence)  for  a  beam  with  a  vertical  distributed  load  proportional  to 
x  +  u  and  an  axial  distributed  load  P.  The  beam  is  allowed  to  be  extensional 
in  that  elongation  of  the  elastic  axis  is  allowed.  The  important  special 
case  where  the  elastic  axis  is  inextensional  is  also  considered.  The  position 
vector  for  a  point  at  a  distance  f;  from  the  elastic  axis  is 

r  =  (x  +  u)I  +  wR  +  £k  (83) 

and  the  axial  strain,  neglecting  elongations  with  respect  to  unity  (ref.  22), 
is 
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Application  of  Hamilton’s  principle  yields 


where  K 
equations 


+  w'6w')  +  EIw" [ (1  +  w’2)6w"  +  w"w'5w'] 


-  P6u  -  mft2(x  +  u)6u  -  K(x  +  u)5wjdx  *  o  (85) 

is  a  given  constant,  A  =  JJ* dA  and  1  =//  £2dA.  The  Euler-Lagrange 
and  boundary  conditions  follow  immediately: 


mft2(x  +  u)  =  P 


[EIw" (1  +  w'2)]"  -  (EIw"2w')'  =  K(x  +  u) 
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EIw" 
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x  =  L 
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(EIw’’)’  =  0  j 


(86a) 

(86b) 

(86c) 


Let  us  consider  for  the  moment  the  case  of  an  inextensional  beam.  To  special¬ 
ize  the  equations  we  must  do  two  things:  (1)  solve  equation  (86a)  for 
EA[u'  +  (w'2/2)]  and  substitute  the  result  in  equation  (86b);  and  then 
(2)  solve  for  u  in  terms  of  w  and  substitute  into  equation  (86b).  Although 
u'  +  (w'2/2)  =  0  for  inextensibility,  EA  -*  00  for  this  case  and  the  product 
EA[u'  +  (w'2/2)]  is  finite,  and  nonzero  in  general.  For  step  (1) 


Ea(u' 
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Step  (2)  yields 


or,  for  EA  ->  °° 
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Thus,  the  single  w  equation  for  an  inextensional  beam  becomes 


with  boundary  conditions  as  above. 

We  now  return  to  the  extensional  case  and  make  use  of  the  quasi-coordinate 
ue  so  that 
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dx 
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Substitution  into  Hamilton's  principle,  equation  (85),  yields 
')  +  EIw" [ (1  +  w'2)6w"  +  w"w'6w']  -  P6u 

-  mfi2|^ 

The  Euler-Lagrange  equations  are 


rL( 

J  <EAUg(Su'  +  w'6w' 

!^x  +  ue  -  -j  J  w'2dx^5u  -  K^x  +  Ug  -  ^  J  w,zdXj^6wl 


dx  =  0  (92) 


-(EAu;)'  -  mft2(x  +  ue  -  f  w’2dx^  =  P 
'  *0  / 

-  (EAu^w' )  '  +  [EIw"  (1  +  w'2)]"  -  (EIw"2w')'  =  +  ug  -  w'2dx^ 


(93a) 


(93b) 


with  boundary  conditions  as  above.  The  underlined  term  in  equations  (92) 
and  (93b)  corresponds  to  (9H/8ug)w' 6w'  in  the  Hamiltonian  (see  equation  (78)). 
Equations  (93a)  and  (93b)  are  identical  to  the  result  of  substituting  equa¬ 
tions  (91)  into  equations  (86a)  and  (86b).  For  the  case  of  inextensibility, 
we  still  must  follow  two  steps:  (1)  solve  equation  (93a)  for  EAu^  and  sub¬ 
stitute  into  equation  (93b);  and  (2)  substitute  ue  =  0  into  equation  (93b). 
The  result  is  equation  (90). 
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To  express  Hamilton's  principle  in  terms  of  ue  and  w  alone  requires 
substitution  of 


6u  =  Su. 


-  f  w'6w'< 
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6u'  =  5ug  -  w'<5w' 


yielding 
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+  [EIw" (1  +  w'2)]" 


The  Euler- Lagrange  equations  are 

-(EAUg)'  -  mfi2^x  +  ue  “  ^  J  w'2dx^  =  P  (96a) 

|  JT  jmfi2^Xj  +  ue  -  j  £  1  w'2dx2^  +  pjdxjj  +  [EIw" (1  +  w’2)]" 

-  (EIw"V)'  =  K^x  +  ue  -  w'2dx^  (96b) 

For  the  case  of  inextensibility,  we  simply  set  ue  =  0  in  equation  (96b)  and 
obtain  equation  (90)  directly. 

Although  no  EA  terms  appear  in  equation  (96b),  it  is  still  suitable  for 
use,  along  with  equation  (96a),  when  EA  is  finite.  The  reason  for  this  is 
now  clear:  when  the  analysis  is  formulated  in  terms  of  the  quasi-coordinate 
ue  at  the  outset,  the  substitution  for  the  tension  EAug  is  automatically 
taken  care  of.  The  rest  of  equation  (96a)  is,  in  effect,  integrated  and  sub¬ 
stituted  automatically.  This  example  explains  why  no  EAug  terms  appear  in 
the  v  and  w  equations  of  references  19-21.  They  do  appear  in  the  v  and  w 
equations  of  reference  8  because  the  case  of  the  tension  axis  being  offset  from 
the  elastic  axis  is  considered. 

At  this  point,  it  seems  appropriate  to  comment  on  the  advantages  of  using 
the  quasi-coordinates  ue  and  4>.  Both  variables  tend  in  some  cases  to  simplify 
the  derivations  of  actual  equations  of  motion  (refs.  3-5,  8,  and  19-21).  It 
is  also  evident  that  taking  the  limit  for  infinite  axial  or  torsional  rigidity 
may  be  somewhat  simpler. 
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6.  IN  THE  LIMIT  OF  INFINITE  TORSIONAL  RIGIDITY 
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Another  way  of  looking  at  the  question  of  the  equivalence  of  the  equations 
in  references  3-5  and  the  two  sets  of  equations  in  reference  8,  demonstrated  in 
section  3,  is  to  consider  the  special  limiting  case  of  infinite  torsional 
rigidity.  It  is  clear  that  no  matter  what  variables  are  used  in  the  analysis 
for  torsion  there  is  a  unique  set  of  equations  governing  bending  deflections 
v  and  w  when  GJ  tends  to  infinity.  For  that  case,  all  GJ  terms  must 
vanish  from  the  v  and  w  equations  in  order  to  preserve  the  proper  mathemati¬ 
cal  structure  of  the  equations.  This  is  similar  to  the  result  of  the  last 
section  in  which  all  EA  terms  vanished  for  the  special  case  of  infinite  EA. 

For  equations  written  in  terms  of  <)>,  such  as  those  of  references  3-5, 
there  are  no  GJ<J> *  terms  in  the  v  and  w  equations  just  as  there  are  no 
EAug  terms  in  the  v  and  w  equations  of  references  19-21.  With  the  order¬ 
ing  scheme  of  references  3-5,  all  of  the  terms  from  the  f  equation  that  are 
integrated  and  substituted  into  the  v  and  w  equations  automatically,  through 
the  process  of  using  the  quasi-coordinate  $,  turn  out  to  be  negligible  (i.e., 
as  stated  in  reference  3,  the  torsion  moment  is  one  order  of  e  smaller  than 
the  bending  moments).  For  infinite  GJ,  the  angle  of  twist  per  unit  length  is 
simply  the  pretwist  component  0£;  the  elastic  twist  <j>  and  its  variation  6d> 
are  simply  set  equal  to  zero  because  elastic  twist  about  the  elastic  axis 
cannot  occur. 


For  equations  written  in  terms  of  0j  or  02,  the  infinite  torsional 
rigidity  constraint  is  more  involved.  If  0j  or  02  were  simply  set  to  zero 
nonzero  GJ  coefficients  would  remain  and  the  equations  would  break  down. 
Instead,  the  GJ  coefficient,  equal  to  -  0*,  must  be  set  equal  to  zero 
and  used  to  eliminate  0.  or  02.  The  following  kinematic  relations  hold  to 
second  order  (see  eqs.  (46)): 


lim  0 1 
GJ-*» 


lim  02 
GJ-**1 


x 

v'Vdx  +  0(e3) 


v'w"dx  +  0(e3) 


I 


(97) 


Then,  the  elastic  torsion  moment  GJ(k^  -  0t)  which  is  finite  and  not  neces¬ 
sarily  zero,  must  be  determined  from  the  torsion  equation  and  substituted  into 
the  v  and  w  equations.  In  this  way  all  or  02  terms  and  all  GJ  terms 
are  eliminated  from  the  v  and  w  equations.  Finally,  if  the  aerodynamic 
pitching  moment  is  not  neglected,  the  virtual  rotation  in  terms  of  v  and  w 
alone  must  be  determined  from  equations  (24)  or  (38)  and  (97): 


J»x 

(w"6v'  -  v"6w')dx  +  0(e3)  (98) 

0 
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for  either  0j  or  02.  The  total  effect  of  pitching  moment  is  then  in  the  terms 
-(w"/«x  dxj)'  in  the  v  equation  and  (v"JLMx  dXj) '  in  the  w  equation. 

It  may  be  shown  that  these  operations  will  yield  a  single  set  of  equations  in 
v  and  w  regardless  of  which  variable  0j  or  02  is  used.  The  only  second- 
order  term  is  the  integral  term  correction  to  the  pitch  angle,  the  "kinematic 
pitch  rotation"  term,  first  identified  in  reference  3.  The  other  places  that 
the  integrals  appear  are  in  third-  and  higher-order  terms  and  not  particularly 
significant.  Therefore,  it  is  not  correct  to  simply  set  0j  or  02  equal  to 
zero  for  general  nonlinear  analysis  of  torsionally  rigid  beams.  In  fact,  when 
any  physical  displacement  variable  u,  v,  w,  0j,  02  that  is  part  of  nonlinear 
Euler-Bernoulli  beam  equations  must  be  eliminated  because  of  infinite  axial, 
bending  or  torsion  stiffness,  it  must  be  done  by  a  process  similar  to  the  one 
outlined  above.  We  have  shown  that  u  4  0  when  EA  -*•  °°,  and  0j  t  02  ^  0 
when  GJ  ■>  °°.  It  is  just  as  easily  shown  that  neither  v  nor  w  vanishes  when 
the  bending  stiffness  in  that  direction  tends  to  infinity.  A  single  set  of 
equations,  consistent  to  second  order,  for  a  torsionally  rigid  rotating  beam 
can  be  obtained  by  simply  setting  tf>  =  0  in  references  4  and  5  or  by  going 
through  the  process  outlined  above  for  either  set  of  the  equations  in  refer¬ 
ence  8. 

Another  result  of  equations  (97)  is  that  three  angles  are  required  to 
describe  the  general  orientation  of  a  reference  frame  in  space.  Thus,  even 
for  infinite  torsional  rigidity  it  is  possible  for  beam  bending  deformations 
to  place  beam  elements  in  such  a  position  that  transformations  based  on  only 
two  angles  cannot  describe  the  orientation  of  a  beam  element.  This  is  seen  in 
the  example  problem  of  section  4,  in  which  three  nonzero  angles  may  be  required 
to  describe  the  orientation  of  a  beam  element,  despite  the  fact  that  *  0. 
The  required  third  angles  are  given  to  second  order  by  equations  (97). 


7.  RIGID  BEAM  APPROXIMATIONS 


When  cantilevered  elastic  beams  are  modeled  by  an  approximate  system 
that  substitutes  a  rigid,  hinged,  spring-restrained  rod  for  the  elastic  beam, 
three  degrees  of  freedom  are  often  considered:  C,  an  angular  deflection  in 
the  x-y  plane  (when  0  =  0  =  0)  or  lead-lag;  0,  an  angular  deflection  out 
of  the  x-y  plane  (when  C  =  0  =  0)  or  flap;  and  0,  the  pitch  angle  (when 
z,  =  8  =  0).  The  nonlinear  equations  for  this  case  are  especially  interesting 
because  six  different  hinge  arrangements  can  be  defined  which  will  lead  to  six 
different  sets  of  nonlinear  equations  (the  linear  equations  being  the  same). 
Simpler  models,  assuming  no  pitch  degree  of  freedom,  have  been  used  in  pre¬ 
liminary  investigations  of  helicopter  rotor  blade  stability  (e.g.,  ref.  23). 
For  this  case  there  are  two  possibilities:  flap-lag,  where  the  lead-lag  hinge 
flaps  with  the  blade,  and  lag-flap,  where  the  flap  hinge  lead-lags  with  the 
blade.  The  two  hinge  arrangements  define  two  different  physical  systems  with 
different  dynamic  properties  although  the  differences  are  not  large.  Differ¬ 
ences  in  the  stability  of  these  two  systems  were  identified  and  discussed  in 
reference  3,  and  reference  6  added  further  numerical  results  and  extensive 
discussion.  It  was  concluded  in  reference  6  that  two  different  sets  of 
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equations  for  a  single  elastia  beam  would  result  from  different  rotation 
sequences  and  that  this  could  be  expected  to  influence  the  stability  of  the 
elastic-beam  model.  This  incorrect  expectation  was  based  on  the  correct  obser¬ 
vation  that  rigid-beam  models  with  two  different  hinge  arrangements  exhibited 
differences  in  stability.  In  fact,  however,  a  single  set  of  equations  for  an 
elastic  beam  can  be  used  to  develop  the  exact  equations  for  a  rigid-beam  model 
with  any  hinge  arrangement.  In  so  doing  one  must  carefully  consider  the 
dependence  of  the  v  and  w  displacements  on  the  hinge  arrangement.  Further¬ 
more  it  can  be  shown  that  the  sequence  of  rotations  used  to  develop  the  equa¬ 
tions  is  arbitrary  for  any  rigid-beam  hinge  arrangement  as  well.  For  simplic¬ 
ity  we  will  perform  these  operations  for  a  second-order  analysis. 

Let  us  assume  first  a  lag-flap  hinge  arrangement  with  the  lead-lag  hinge 
at  the  root  and  the  flap  hinge  offset  from  the  root  at  axial  distance  e. 

The  limit  must  be  taken  as  e  ->  0  to  recover  the  correct  hinge  arrangement 
even  for  the  coincident  hinge  model  used  in  reference  6.  The  following  dis¬ 
placement  functions  are  assumed 


v  =  £x  ,  v'  =  ;  ,  v"  =  0  x  >  0 


(99a) 


but,  so  that  the  flap  hinge  lead-lags  with  the  beam. 


\ 


x  >  0 


> 


(99b) 


w"  =  lim  3<5 (e)  =86(0)  x  >  0 
e-*0 


where  6(e)  and  6(0)  refer  to  the  Dirac  delta  function.  For  a  flap-lag 
arrangement,  we  need  to  reverse  the  expressions  in  equations  (99)  so  that 


w  *  8x  ; 

w’  =  6  ; 

V  =  £X 

x  >  0 

1° 

x  =  0 

v- 1 

x  >  0 

v"  «  56  (0) 

x  >  0 

x  >  0 


(100) 


To  reduce  the  elastic-beam  equations  in  references  4  and  5  to  the  torsionally 
rigid  case  we  simply  set  $  ■  0.  The  resulti  equations  would  then  contain 
the  integrals  in  equations  (97)  that  have  different  values  for  the  two  hinge 
arrangements  defined  by  equations  (99)  and  (100).  Evaluation  of  the  integrals 
for  the  lag- flap  and  flap- lag  hinge  arrangements  yields 
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atomaKU' 


for  lag-flap: 


for  flap-lag: 


fX  v’V 

^  n 


dx  =  0 


(101) 


f  v*w"dx  =  r  66(0)5  dx  =  65 

'0  *0 

f  v’Vdx  =  f  56(0)6  dx  =  65 

*0  '0 

Jv'w"dx  *  0 

n 


(102) 


The  differences  between  equations  (101)  and  (102)  result  in  two  different  sets 
of  equations  for  the  two  different  hinge  arrangements,  and  these  equations 
are,  to  second  order,  identical  to  the  rigid-blade  equations  derived  in 
reference  6. 

For  two  degrees  of  freedom  6  and  5  and  a  given  hinge  arrangement,  any 
sequence  of  rotations  (e.g.,  5j,  6j,  0j  or  62*  £2*  e2^  can  used  to  formu¬ 
late  the  equations  although  it  would  certainly  be  more  natural  to  use  the 
angles  associated  with  rotations  about  the  hinges.  The  reconciliation  of 
three  angles  with  only  two  degrees  of  freedom  is  made  by  use  of  a  holonomic 
constraint  equation  that  relates  0  to  6  and  5  (ref.  10,  p.  55).  For 
example,  consider  writing  equations  for  a  flap- lag  hinge  arrangement  in  terms 
of  5lf  6j,  0,,  instead  of  in  terms  62»  52  (0?  not  needed  because  the 
axes  about  which  02  and  52  occur  are  along  the  hinge  axes).  Thus  02  =  0 
and  equation  (40)  yields 


02  =  0 


+  tan 


"( tt) 


(103) 


0j  *  -tan-1(sin  6j  tan  5j)  =  -6^  +  0(e3) 


(104) 


In  fact,  as  long  as  the  proper  constraint  is  accounted  for,  any  set  of  vari¬ 
ables  can  be  used  to  write  the  equations  for  any  physical  hinge  arrangement. 
When  the  constraints  are  properly  accounted  for,  the  equations  of  motion  for  a 
given  physical  hinge  arrangement  are  equivalent,  although  different  in  appear¬ 
ance,  for  any  sequence  of  rotations  used  to  describe  the  beam  orientation. 

The  dynamic  behavior  of  the  system  is  then  independent  of  the  transformation 
sequence  used  in  describing  the  orientation  of  the  blade.  The  equations  and 
dynamic  behavior  do,  however,  depend  on  the  physical  arrangement  of  the  hinges. 
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8.  CONCLUSIONS 


In  this  report  the  large  deformation  geometry  for  an  Euler-Bernoulll  beam 
has  been  developed.  Some  aspects  of  the  beam  kinematics  discussed  in  refer¬ 
ences  3-8,  12,  16,  and  17  have  been  clarified.  The  following  points  summarize 
the  results  of  this  study: 

1.  The  large-deformation  kinematics  for  an  Euler-Bernoulli  beam  are 
developed,  including  the  transformation  matrix  relating  the  local  principal 
axes  in  the  deformed  state  to  space-fixed  Cartesian  axes,  the  components  of 
angular  velocity  and  virtual  rotation  vectors,  the  torsion,  and  the  components 
of  bending  curvature.  The  values  of  all  the  geometric  quantities  are  unique 
at  a  given  instant  in  time,  but  the  form  of  the  expressions  themselves  may 
depend  on  the  variables  used  in  describing  the  orientation  of  the  cross  sec¬ 
tion  during  deformation. 

2.  The  angles  5j,  Bj,  0j  are  not  the  same  as  the  angles  82,  C2,  ®2» 
both  sets  of  which  are  discussed  in  the  text.  Although  the  geometry  is  unique, 
the  use  of  the  different  sequences  of  rotations  may  change  the  appearance  of 
the  geometric  expressions. 

3.  The  exact  expressions  for  the  large-deformation  geometry  that  are 
derived  in  this  report  do  not  depend  on  the  sequence  of  rotations  used  in 
defining  the  beam  cross-section  orientation.  Some  of  these  quantities  were 
also  derived  previously  in  references  3-5. 

4.  The  stability  and  dynamic  behavior  of  an  elastic  beam  are  not  depen¬ 
dent  on  the  choice  of  angles  used  to  describe  the  orientation  of  the  local 
beam  cross  section  during  deformation.  In  fact,  the  two  sets  of  equations 
derived  in  reference  8,  based  on  two  different  sets  of  angles,  are  shown  to  be 
equivalent  to  each  other  and  to  the  equations  of  references  4  and  5. 

5.  The  exact  relationship  of  the  pitch  rotation  angles  8.  and  ©2  and 
the  elastic  twist  angle  <f>,  established  in  reference  3,  is  verified  along  with 
the  exact  transformation  matrix  first  derived  therein.  The  integrals  asso¬ 
ciated  with  these  quantities  are  correct  as  originally  derived  in  reference  3, 
and  serve  to  show  that  there  can  be  a  kinematic  pitch  rotation  of  the  beam 
cross  section  that  depends  on  bending  alone,  even  when  there  is  no  pretwist 
and  no  elastic  twist. 

6.  The  kinematic  pitch  rotation  due  to  bending  only  that  may  occur  even 
in  torsionally  rigid  beams  is  seen  to  be  mathematically  analogous  to  a  kine¬ 
matic  shortening  that  may  occur  even  in  axially  rigid  beams  due  to  bending 
only.  The  torsion  variable  <f>  used  in  this  report  is  the  elastic  angle  of 
twist,  which  is  devoid  of  kinematic  pitch  rotation  due  to  bending  alone.  Its 
mathematical  analog  is  the  axial  deflection  due  to  longitudinal  strain  of  the 
elastic  axis  ue,  which  is  devoid  of  kinematic  axial  deflection  (foreshorten¬ 
ing)  due  to  bending  only. 

7.  Both  the  elastic  twist  angle  <p  and  the  axial  deflection  due  to 
longitudinal  strain  of  the  elastic  axis  ue,  because  they  are  related  to 


physical  deflections  through  integrals  and  have  integrals  of  angular  velocity 
and  velocity  components  for  kinetic  analogs,  belong  to  a  certain  class  of 
variables  called  quasi-coordinates.  They  have  both  been  used  in  Hamilton's 
principle  without  the  rigorous  justification  provided  in  this  report.  Although 
not  previously  established,  their  use  in  references  4,  5,  8,  and  19-21  is 
found  to  be  valid. 

8.  Three  angles  are  shown  to  be  required  in  general  to  describe  the  local 
cross-section  orientation  for  an  elastic  beam,  even  if  the  torsion  stiffness 

is  infinite.  Only  one  set  of  equations  exist  for  bending  deflections  v  and 
w  of  an  inextensional,  torsionally  rigid  Euler-Bernoulli  beam,  regardless  of 
the  three  angles  used  in  the  derivation. 

9.  When  the  proper  mathematical  limits  for  different  hinge  arrangements 
in  a  hinged  rigid-beam  model  are  imposed  on  the  elastic-beam  equations  with 
rigid-body  mode  shapes,  the  correct  hinged,  rigid-beam  equations,  for  each 
particular  hinge  arrangement,  are  obtained. 

10.  Any  three  angles  can  be  used  to  describe  the  orientation  of  any 
hinged,  rigid  body  as  long  as  the  constraints  that  relate  the  rotations  about 
the  hinges  to  the  variables  in  the  derivation  are  properly  accounted  for.  Of 
course,  it  may  be  more  natural  to  use  the  physical  angles  as  variables.  It  is, 
therefore,  not  the  transformation  sequence  that  influences  the  dynamic 
behavior  —  it  is  the  physical  arrangement  of  the  hinges. 

11.  The  bending  curvatures,  torsion,  and  [T]  matrix  are  developed  in  the 
appendix  for  six  possible  sequences  of  the  angles  ?,  8,  and  0.  Since  values 
of  the  torsion  and  elements  of  the  [T]  matrix  are  fixed  at  a  given  Instant  in 
time,  a  change  of  variable  may  be  found  from  either  quantity  to  relate  the 
geometry  from  any  one  sequence  to  that  of  any  other. 
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APPENDIX  A 


KINEMATICS  FOR  OTHER  SEQUENCES 


In  the  text,  the  development  was  restricted  to  only  two  of  the  six  pos¬ 
sible  sequences  of  rotations  C»  6,  0.  In  the  appendix  we  will  present  expres¬ 
sions  for  the  [T]  matrix,  the  bending  curvatures,  and  the  torsion  for  all  six 
possible  sequences  of  rotations:  (1)  Cj,  0j,  Q  j ;  (2)  02,  ^2»  02»  (3)  C3,  03, 
e3;  (4)  04»  e4,  c4;  (5)  e5,  C5.  03;  and  (6)  06,  B6»  C6>  The  necessary  changes 
of  variable  to  demonstrate  equivalence  can  be  formulated  from  the  torsion  and 
then  the  elements  of  [T]  can  be  shown  to  be  equivalent  for  all  six  sequences. 
Alternatively,  the  change  of  variable  can  be  formulated  from  the  T23  element 
of  [T]  and  then  the  torsion  and  bending  curvatures  can  be  shown  to  be  equiva¬ 
lent  for  all  six  sequences.  It  is  the  second  approach  that  is  employed  here. 

We  will  make  use  of  equations  (5)  extensively,  which  we  repeat  here  for 
convenience 


(x  +  u)+  =  Tn  \ 

v+  =  T12  i  (Al) 

w+  -  T13  J 


These  relationships  hold  for  all  possible  sequences  and  types  of  rotations. 
The  [T]  matrix,  the  bending  curvatures,  and  the  torsion  are  now  presented  for 
each  of  the  six  sequences. 


(1)  Cj,  Bj,  0J--  The  [T]  matrix  is 


[T] 


c  c 
8l  «1 


"^Al  '  S01S31CCi  C©Al  ‘  “cAAl 
"c0xsBiCCi  +  8eAl  'S6iC?i  '  S£iSBiC0i 


The  bending  slopes  are,  from  equation  (Al) 


v+  =  c0  s  ;  v'  -  Cl  +  0(e3) 
Cl 


w+  ■  s 


Si 


;  w'  *  Bi  +  0(e3) 


Si 


C.  fi 

c»i 


c6ice: 


(A2) 


(A3) 


The  rotation  rate  vector  is 


K  -  Cl*  -  st(cr,  J  -  s5]i)  +  e+i  - 


(A4) 
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+  <jJ  +  <kk 


Thus 


tc.  =  et  +  ?+s  *  8!  +  v"w'  +  0(e3) 

1  -i  Pi  1 

Kj  *  "6tcei  +  ?tcBlsei  “  “W”  +  91V"  +  0(e3) 

Kk  =  6ts6l  +  4%%  =  V"  +  91W"  +  °(e3)  J 

(2)  $2>  ^2»  02*“  T^e  is 


[T]  = 


cr  cft 
^2  02 


?2 


cr  SR 
^2  B2 


(A5) 


'Ce2C02S^2  S^2S02  C02C?2  **2*2  " 


‘2  h 2  M2 


"Se2C02  +  S92S?2CB2 


"S02CC2  S02Se2S?2  +  CB2C92 


The  bending  slopes  are,  from  equations  (Al) 

;  v'  =  ?2  +  0(e3) 


v+  =  s. 


w+  =  cr  s«  5 
52  B2 


w'  =  &2  +  0(e  3) 


(A6) 


(A7) 


The  rotation  rate  vector  is 


‘  '  -e"J  +  5t(V  ■  V)  +  e*‘  '  Ki1  +  ‘j1  +  \k 


(A8) 


Thus 


i  =  e2  -  0*s  -  0'  -  v'w"  +  0(e3) 

4  *  -etcQ  c.  +  ^2se  55  ~w"  +  e?v"  +  0(e3)! 


2  02  ’2 


*k  =  4C92  +  Bt802^2  "  V"  +  02w"  +  °<83> 


(A9) 
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•aT>  ’-  ^ 


(3)  c3,  03,  8 3« -  The  [T]  matrix  is 


[T] 


c83%  +  S83S03%  ce3%  "  se3se3%  ce,se 


"cfl  sr 
9  3  ^  3 


cft  cr 
9  3  ^  3 


-C_  S .  +  C -  s„  s.  -c„  s„  c  -  s 


c3  e3  8,  e,  ? 


3  3  ’  3 


63  03  ?3  6,  C 


3  K3 


Vce. 


3  s3  3  v3 


(A10) 


The  bending  slopes  are,  from  equation  (Al) 

V+  =  cr  Sr  "  s o  sfl  cr  » 

93  ^  3  93  9  3  ’ 3 


v'  -  c3  -  63e3  +  0(e3) 


W+  =  c„  s 


9 3  93 


w'  =  e3  +  0(e3) 


(All) 


The  rotation  rate  vector  is 

K  = 


+  *  \*)  -  6tJ  *  'i1  +  ‘ji  +  'kk 


(A12) 


Thus 


Ki  =  4c03sb3  +  etce3  =  e3  +  v"w'  +0(e3) 

Kj  -  -6+  +  ?3s0  *  -w"  +  03v"  +  0(e3) 

Kk  "  4cb3c03  -  9ts83  =  V"  +  9 3W"  +  °(e3) 


(A13) 


(4)  64.  04»  C4*-  The  [T]  matrix  is 


[T] 


cr  CR  “  Sr  Sfi  SR 
S  ^  *r»  ^  P  P 


Sr  Cfi  Cr  SR  +  sr  Sfl  CR 

**4  94  ^4  B4  ^4  94  94 


S.  c„  -  s„  c  sn  c,  cQ  -s.  s^  +  c„  s.  c 


c4  e4  e4  c4  e4 


"V^ 


^4  9 4 


-Sr 


c4  c4  e4  e4 

W 


(A14) 


The  bending  slopes  are,  from  (Al) 

v+  *  ar  CR 
C4  e4 


w+  =  C  s  +  s  s.  c  ;  w'  *  8t 
44  b4  44  °4  p4 


+  0(e3)  | 

+  C404  +  0(e3)  I 


(A15) 
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The  rotation  rate  vector  is 


‘  -  -<•>  +  *  -,/)  *  <>  ■  «!*■  +  «jJ  * 


Thus 


K  -  etc  -  ets.  c  -  e'  -  vV  +  o(c3) 

1  *4  ^4  4 

kj  "  -etc  c  -  eta  *  -w"  +  e4v"  +  0(e3) 

J  ^4  u4 


Kk  “  +  6tse  =  v"  +  e4w"  +  0(e3) 

a  4 


(A16) 


(A17) 


(5)  05,  ?5,  6-.-  The  [T]  matrix  is 


[T] 


Cft  Cr  CR  S r  C A  ■  SR  SA  CA  Sr  SA  +  S«  CA 

e5  ?5  35  cs  e5  b5  e5  b5  c5  e5  Be  e 


5  5 


cr  CA 
?5  05 


Cr  s0 
C5  05 


-%%  -  \\\  \\  -  \\\ 


The  bending  slopes  are,  from  equations  (Al) 


v+  =  cQ  s„  c0  -  sQ  sQ  -  C5  -  6565  +  0(e3) 


'*5  ^ 5  9 5  65  6 5 


w+  =  CA  Sr  SA  +  Sft  CA  =  S5  +  C505  +  0(e3) 
6S  «5  05  SS  95 


(A18) 


(A19) 


The  rotation  rate  vector  is 

5-  ’5Vce 

Thus 


'  ‘  9-r  +  «t(ce,‘  -  YJ'  - 


)  -  8tJ  •  “i1  +  "l1  + 


k.  =  etcfi  c  +  ctsa  *  9 1  +  v’V  +  0(e 3) 

1  5  ?5  5  B5  5 

k.  -  -et  -  etsr  =  -w"  +  e.v"  +  o(e3) 

j  5  5  s5  5 

"k  "  4ce,  "  et%V  "  v"  +  95w"  +  0(e3) 

D  5 


(A20) 


(A21) 
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(6)  0g,  gg,  Cg.-  The  [T]  matrix  is 


[T] 


C6  96  ?6  9 6  C6  96  9 6  C6  ®6  ?6  96  06 

-s  c  c  C -  +  s  s  s.  c  s.  -  s  s  c 

9 6  66  5 6  \  ®6  C6  06  ?6  *6  86 


-s. 


V'e 


(A22) 


The  bending  slopes  are,  from  equation  (Al) 

v+  ■  SC.V  -  C£CVV  !  v’  •  s6  -  S6»6  +  0<s3) 

ob  6  6  6 

w+  =  sr  Se  +  Cc  SB  CA  ’  w’  =  e6  +  ^6e6  +  0(e3) 
^6  6  ’6  p 6  °6 


(A23) 


The  rotation  rate  vector  is 


=  etx  *  9t  (ce6J  +  se6K)  +  ^  =  KiA  +  <ji  +  Kkk 


(A24) 


Thus 


Ki  =  6tv&g  -  =  9 6  -  vV'  +  (e3) 

Kj  *  “6t%  *  6tS5gCBg  “  “W"  +  e6V"  +  °(e3) 

Kk  =  4  "  96S6c  *  V"  +  96W”  +  °(e3) 

b 


(A25) 


Demonstration  of  Equivalence 


It  can  be  shown  that  changes  of  variables  constructed  from  equating  three 
elements  of  the  [T]  matrix  for  any  two  sequences  will  lead  to  identical  expres¬ 
sions  for  bending  curvatures  and  the  torsion.  Only  one  example  is  presented 
for  illustration  —  the  equivalence  of  (1)  and  (2).  For  convenience  we  choose 
to  obtain  the  change  of  variable  from  T23.  From  equations  (Al),  (A2), 
and  (A6) ,  the  following  relations  hold: 


v+  =  CR  Sr  "  3r 
91  C1  C2 


WT  *  s  =  c  s 

61  C2  S2 


T2  3  =  CR  SA  =  SA  CA  '  C  A  Sr  SR 
3  9  1  91  9 2  9 2  9 2  C2  6 2 


(A26) 
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Substitution  of  equations  (A28)  into  (A27)  yields 


ft+  _  ?2%SB2 

t  ‘K‘K 

/+  +  ^ 

+  ?+  %%  pCS2C62  '  ^lS^2S^2 

62  32  \  1~ClZS&2  > 


B2Sr  +  ?2Sfi  cr  CP 
a+  ft+  g2  g2  g2  p2 

61  °2  "  .  2  2 

1  -  s„ 

g2 

Substitution  of  equations  (A27)  through  (A29)  into  equations  (A5)  yields  for 


s+.  +  &„  <=  c  %ck,  . 

K.  =  9,  -  -  +  - - -  S, 

*>2 

£  0  P-3 


-  e+  -  s  4 
^2 


52  +  sr  sr 

2  e2  c2 


o+  +  > 

32Cs2C62  C2S^2SB2 
1  -  C?  s£ 


^2  g2 


(A29) 


which  is  identical  to  the  expression  for  in  equations  (A9).  The  bending 

curvatures  also  transform  identically.  Similarly,  the  other  four 
sequences  (3)-(6)  can  be  shown  to  be  equivalent  to  (1).  The  changes  of 
variable  to  second  order  may  be  determined  by  inspection,  which  proves  that 
all  six  sequences  are  also  equivalent  to  second  order. 
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